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GENERAL POTENTIAL THEORY OF ARBITRARY WING SECTIONS 
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SUMMARY 

This report gives an exact treatment of the problem of 
determining the 2-dimensional potential flow around 
wing sections of any shape. The treatment is based 
directly on the solution of this problem as advanced by 
Theodorsen in N. A. C. A. Technical Report No. Ul. 
The problem condenses into the compact form of an inte- 
gral equation capable of yielding numerical solutions by 
a direct process. 

An attempt has been made to analyze and coordinate 
the results of earlier studies relating to properties of wing 
sections. The existing approximate theory of thin wing 
sections and the Joukowsky theory with its numerous 
generalizations are reduced to special cases of the general 
theory of arbitrary sections, permitting a clearer perspec- 
tive of the entire field. The method not only permits the 
determination of the velocity at any point of an arbitrary 
section and the associated lift and moments, but furnishes 
also a scheme for, developing new shapes of preassigned 
aerodynamical properties. The theory applies also to 
bodies that are not airfoils, and is of importance in other 
branches of physics involving potential theory. 

INTRODUCTION 

The solution of the problem of determining the 
2-dimensional potential flow of a nonviscous incom- 
pressible fluid around bodies of arbitrary shape can be 
made to depend on a theorem in conformal represen- 
tation stated by Riemann almost a century ago, 
known as the fundamental theorem of conformal rep- 
resentation. This theorem is equivalent to the state- 
ment that it is possible to transform the region 
bounded by a simple curve into the region bounded by 
a circle in such a way that all equipotential lines and 
stream lines of the first region transform respectively 
into those of the circle. The theorem will be stated 
more precisely in the body of this report and its sig- 
nificance for wing section theory shown — suffic 5 it at 
present to state that if the analytic transformation by 
which the one region is transformed conformally into 
the region bounded by the circle is known, the poten- 
tial field of this region is readily obtained in terms of 
the potential field of the circle. 

A number of transformations have been found by 
means of which it is possible to transform a circle into 



a contour resembling an airfoil shape. It is obviously 
true that such theoretical airfoils possess no particular 
qualities which make them superior to the types of more 
empirical origin. It was probably primarily because 
of the difficulty encountered in the inverse problem, 
viz, the problem of transforming an airfoil into a 
circle (which we shall denote as the direct process) 
that such artificial types came into existence. The 
2-dimensional theoretical velocity distribution, or what 
is called the flow pattern, is known only for some 
special symmetrical bodies and for the particular class 
of Joukowsky airfoils and their extensions, the out- 
standing investigators 1 being Kutta, Joukowsky, and 
von Mises. Although useful in the development of 
airfoil theory these theoretical airfoils are based solely 
on special transformations employing only a small 
part of the freedom permitted in the general case. 
However, they still form the subject of numerous 
isolated investigations. 

The direct process has been used in the theory of 
thin airfoils with some success. An approximate 
theory of thin wing sections applicable only to the 
mean camber line has been developed 2 by Munk and 
Birnbaum, and extended by others. However, at- 
tempts 3 which have been made to solve the general 
case of an arbitrary airfoil by direct processes have 
resulted in intricate and practically unmanageable 
solutions. Lamb, in his "Hydrodynamics" (reference 
I, p. 77), referring to this problem as dependent upon 
the determination of the complex coefficients of a 
conformal transformation, states: "The difficulty, 
however, of determining these coefficients so as to 
satisfy given boundary conditions is now so great as 
to render this method of very limited application. 
Indeed, the determination of the irrotational motion of 
a liquid subject to given boundary conditions is a 
problem whose exact solution can be effected by direct 
processes in only a very few cases. Most of the cases 
for which we know the solution have been obtained by 
an inverse process; viz, instead of trying to find a 
value of <t> or ^ which satisfies [the Laplacian] = 0 
or vV=0 and given boundary conditions, we take 
some known solution of these differential equations 

i See bibliography given in reference 9, pp. 24, 84, and 583. 

i Cf. footnote 1. 

< See Appendix II of this paper. 
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and inquire what boundary conditions it can be made 
to satisfy." 

In a report (reference 2) recently published by the 
National Advisory Committee for Aeronautics a gen- 
eral solution employing a direct method was briefly- 
given. It was shown that the problem could be stated 
in a condensed form as an integral equation and also 
that it was possible to effect the practical solution of 
this equation for the case of any given airfoil.- A 
formula giving the velocity at any point of the surface 
t>f an arbitrary airfoil was developed. The first part 
of the present paper includes the essential develop- 
ments of reference 2 and is devoted to a more com- 
plete and precise treatment of the method, in particu- 
lar with respect to the evaluation of the integral 
equation. 

In a later part of this paper, a geometric treatment 
of arbitrary airfoils, coordinating the results of earlier 
investigations, is given. Special airfoil types have 
also been studied on the basis of the general method 
and their relations to arbitrary airfoils have been 
analyzed. The solution of the inverse problem of 
creating airfoils of special types, in particular, types of 
specified aerodynamical properties, is indicated. 

It is hoped that this paper will serve as a step 
toward the unification and ultimate simplification of 
the theory of the airfoil. 

TRANSFORMATION OF AN ARBITRARY AIRFOIL INTO 
A CIRCLE 

Statement of the problem.— The problem which this 
report proposes to treat may be formulated as follows. 
Given an arbitrary airfoil 4 inclined at a specified angle 
in a nonviscous incompressible fluid and translated 
with uniform velocity V. ' To determine the theoreti- 
cal 2-dimensional velocity and pressure distribution at 
all points of the surface for all orientations, and to 
investigate the properties of the field of flow surround- 
ing the airfoil. Also, to determine the important 
aerodynamical parameters of the airfoil. Of further 
interest, too, is the problem of finding shapes with 
given aerodynamical properties. 

Principles of the theory of fluid flow. — We shall 
first briefly recall the known basic principles of the 
theory of the irrotational flow of a Motionless incom- 
pressible fluid in two dimensions. A flow is termed 
"2-dimensional" when the motion is the same in all 
planes parallel to a definite one, say xy. In this case 
the linear velocity components u and v of a fluid 
element are functions of x, y, and i only. 
' The differential equation of the lines of flow in this 
case is 

v da; — u dy = 0 

* By an airfoil shape, or wing section, is roughly meant an elongated smooth shape 
rounded at the leading edge and ending in a sharp edge at the rear. All practical 
airfoils are characterized by a lack of abrupt change of curvature except for a rounded 
nose and a small radius of curvature at the tail. 



11 = 



and the equation of continuity is 

bx by~ or bx~ by 

which shows that the above first equation is an exact 
differential. 
If Q=c is the integral, then 

by bx 

This function Q is called the stream function, and 
the lines of flow, or streamlines, are given by the equa- 
tion Q = c, where e is in general an arbitrary function 
of time. 

Furthermore, we note that the existence of the 
stream function does not depend on whether the motion 
is irrotational or rotational. When rotational its 
vorticity is 

^ bx by bx 2 by 2 

which is twice the mean angular velocity or "rotation" 
of the fluid element. Hence, in irrotational flow the 
stream function has to satisfy 



dx 2 + d^ -u 



(2') 



Then there exists a velocity potential P and we have 



bP 

dx 

bP 



bQ 
by 



— v=- 



(1) 



(2) 



bQ 

bx 

The equation of continuity is now 

b 2 P,b 2 P_. 
btf + by*' 0 

Equations (1) show that 

dPbQ + dPbQ = 0 
bx bx by by 

so that the family of curves 

P = constant, Q — constant 

cut orthogonally at all their points of intersection. 

For steady flows, that is, flows that do not vary 
with time, the paths of the particles coincide with the 
streamlines so that no fluid passes normal to them. 
The Bernoulli formula then holds and the total pres- 
sure head H along a streamline is a constant, that is 

%p^+p'=H 

where p' is the static pressure, v the velocity, and p 
the density. If we denote the undisturbed velocity 
at infinity by V, the quantities p'—p'-o by p, and 
% p V 2 by q, the Bernoulli formula may be expressed as 



(3) 
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The solutions of equations (2) and (2'), infinite in 
number, represent all possible types of irrotational 
motion of a nonviscous incompressible fluid in two 
dimensions. For a given problem there are usually 
certain specified boundary conditions to be satisfied 
which may be sufficient to fix a unique solution or a 
family of solutions. The problem of an airfoil moving 
uniformly at a fixed angle of incidence in a fluid field 
is identical with that of an airfoil fixed in position and' 
fluid streaming uniformly past it. Our problem is 
then to determine the functions P and Q so that the 
velocity at each point of the airfoil profile has a direc- 
tion tangential to the surface (that is, the airfoil con- 
tour is itself a streamline). and so that at infinite dis- 
tance from the airfoil the fluid has a constant velocity 
and direction. 

The introduction of the complex variable, z—x+iy, 
simplifies the problem of determining P and Q. Any 
analytic function w(z) of a complex variable z, that is, 
a function of z possessing a unique derivative in a 

z Plane 




Figure 1.— Conformal property of analytic [unctions 

region of the complex plane, may be separated into its 
real and imaginary parts as w(z)=w(x+iy)=P(x, y) 
+iQ fr, y), determining functions P and Q which may 
represent the velocity potential and stream function of 
a possible fluid motion. Thus, analytic functions of a 
complex variable possess the special property that the 
component functions P and Q satisfy the Cauchy- 
Riemann equations (eq. (1)), and each therefore also 
satisfies the equation of Laplace (eq. (2)). Conversely, 
any function P(x, y) + iQ(.x, y) for" which P and Q 
satisfy relations (1) and (2) may be written as w{x+ 
iy)=£w(z). The essential difficulty of the problem is 
to find the particular function w(z) which satisfies the 
special boundary-flow conditions mentioned above for 
a specified airfoil. 

The method of conformal representation, a geomet^ 
ric application of the complex variable, is well adapted 
to this problem. The fundamental properties of trans- 
formations of this type may be stated as follows: 
Consider a function of a complex variable z=x+iy, 
say w (z) analytic in a given region, such that for each 
value of z, w(z) is uniquely defined. The function 
w(z) may be expressed as w={+ in where £ and y are 



each real functions of x and y. Suppose now in the 
xy complex plane there is traced a simple curve /(z)- 
(Fig. 1.) Each value of z along the curve defines a 
point w in the w plane and/(s) maps into a curve /(w) 
or F(z). Because of the special properties of analytic 
functions of a complex variable, there exist certain 
special relations between /(z) and F(z). 

, The outstanding property of functions of a complex 
variable analytic in a region is the existence of a unique 
derivative at every point of the region. 

Urn 110—10' , 



dw 
dz" 



or 



'z-»z' z—z' 
dw = pe'Tdz 



This relation expresses the fact that any small curve 
22' through the point z is transformed into a small 
curve im>' through the point w by a magnification p 
and a rotation 7; i. e., in Figure 1 the tangent t will 
coincide in direction with T by a rotation 7 = /3 — a. 
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Ficukb 2.— Orthogonal network obtained by a conformal transformation 

This is also true for any other pair of corresponding 
curves through z and w, so that in general, angles 
between corresponding curves are preserved. In par- 
ticular, a curve zz" orthogonal to zs' transforms 
into a curve ww" orthogonal to vm'. 

It has been seen that an analytic function /(s) may 
be written P(x, y) + iQ{x, y) where the curves P= con- 
stant and Q= constant form an orthogonal system. 
If then /(z) is transformed conformally into /(w) 
= P(f, v ) + iQ& v) that is into .flw(z)] = F(z) - B(x, y) 
+ iS(x, y), the curves P{x, y) = constant, Q(x, y) = con- 
stant map into the orthogonal network of curves 
fife y) = constant, S(x, y) = constant. (Fig. 2.) If the 

~ 0 =p is zero at a point w, the trans- 
formation at that point is singular and ceases to be 
conformal. 

We may use the method of conformal transforma- 
tions to find the motion about a complicated boundary 
from that of a simpler boundary. Suppose w{z) is a 
function which corresponds to any definite fluid motion 
in the z plane, for instance, to that around a circle. 
Now if a new variable f is introduced and 2 set equal 



magnification I -j^ 
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to any analytic function of f, say z=/(f), then w{z) 
becomes )] or W{$) representing a new motion in 
the f plane. This new motion is, as has been seen, 
related to that in the z plane in such a way that the 
streamlines of the z plane are transformed by z=?/(f) 
into the streamlines of the f plane. Thus, the con- 
tour into which the circle is transformed represents 
the profile around which the motion W(f ) exists. The 
problem of determining the flow around an airfoil is 
now reduced to finding the proper conformal transfor- 
mation which maps a curve for which the flow is known 
into the airfoil. The existence of such a function was 
first shown by Biemann. 

We shall first formulate the theorem for a simply 
connected region 5 bounded by a closed curve, and 
then show how it is readily applied to the region 
external to the closed curve. The guiding thought 
leading to the theorem is simple. We have seen that 
an analytic function may transform a given closed 
region into another closed region. But suppose we 
are given two separate regions bounded by closed 
curves — does there exist an analytic transformation 
which transforms one region conformally into the 
other? This question is answered by Biemann's 
theorem as follows: 

Riemann's theorem. — The interior T of any simply 
connected region (whose boundary contains more than 
one point, but we shall be concerned only with regions 
having closed boundaries, the boundary curve being 
composed of piecewise differentiable curves [Jordan 
curve], corners at which two tangents exist being per- 
mitted) can be mapped in a one-to-one conformal 
manner on the interior of the unit circle, and the 
analytic 6 function f=/(z) which consummates this 
transformation becomes unique when a given interior 
point Zo of T and a direction through Zo are chosen to 
correspond, respectively, to the center of the circle and 
a given direction through it. By this transformation 
the boundary of T is transformed uniquely and con- 
tinuously into the circumference of the unit circle. 

The unit circle in this theorem is, of course, only a 
convenient normalized region. For suppose the re- 
gions Ti in the f plane and T 2 in the w plane are 
transformed into the unit circle in the z plane by 
f=/(z) and w=F(z), respectively, then Ti is trans- 
formed into T 2 by f = $(w), obtained by eliminating z 
from the two transformation equations. 

In airfoil theory it is in the region external to a closed 
curve that we are interested. Such a region can be 
readily transformed conformally mto the region in- 
ternal to a closed curve by an inversion. Thus, let us 
suppose a point Zo to be within a closed curve B whose 

• A region of the complex plane is simply connected when any closed contour lying 
entirely within the region may be shrunk to a point without passing put of the region. 
Cf. reference 3, p. 367, where a proof of the theorem based on Green's function is 
given. 

« Attention is here directed to the fact that an analytic function is developable at 
a point in a power series convergent in any circle about the point and entirely 
within the region. 



external region is r, and then choose a constant k 
such that for every point z on the boundary of T, 
\z — z 0 |>&. Then the inversion transformation w = 
k 

— — will transform every point in the external region 

Z Zq 

T into a point internal to a closed region r' lying 
entirely within B, the boundary B mapping into the 
boundary of r', the region at infinity into the region 
near Zo. We may now restate Biemann's theorem as 
follows: 

One and only one analytic function f =/(z) exists by 
means of which the region r external to a given curve 
B in the f plane is transformed conformally into the 
region external to a circle C in the z plane (center at 
z—Q) such that the point z= <=° goes into the point 

f = co and also ^p = l at infinity. This function can 

be developed in the external region of C in a uniformly 
convergent series with complex coefficients of the form 

f -m=/(z)-m.=2+!+5+^+ • • • (4). 

by means of which the radius R and also the constant 
m are completely determined. Also, the boundary B 
of T is transformed continuously and uniquely into the 
circumference of C. 

It should be noticed that the transformation (4) is 
a normalized form of a more general series 

f-m=a 0 +a_ l z+J+J+ 

and is obtained from it by a finite translation by the 
vector do and a rotation and expansion of- the entire 
field depending on the coefficient a_!. The condition 
o_i = l is necessary and sufficient for the fields at 
infinity to coincide in magnitude and direction. 

The constants c< of the transformation are functions 
of the shape of the boundary curve alone and our 
problem is, really, to determine the complex coeffi- 
cients defining a given shape. With this in view, we 
proceed first to a convenient intermediate trans- 
formation. 

a 2 

The transformation t=z'+—. — This initial trans- 

z 

formation, although not essential to a purely mathe- 
matical solution, is nevertheless very useful and 
important, as will be seen. It represents also the key 
transformation leading to Joukowsky airfoils, and is 
the basis of nearly all approximate theories. 

Let us define the points in the f plane by f = x + iy 
using rectangular coordinates (x, y), and the points in 
the z' plane by z' = ae* +i6 using polar coordinates 
(ae*, 0). The constant a may conveniently be con- 
sidered unity and is added to preserve dimensions. 
We have 

r-s'+jP (5) 
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(6) 



and substituting z' = ae# + *' 

we obtain f = 2a cosh (4>+id) 

or £ = 2a cosh if/ cos 0+2ia sinh ^ sin 0 

Since £=z+iy, the coordinates (2, y) are given by 

x= 2a. cosh ^ cos 0 

2/ = 2a sinh ^ sin 0 

If i^=0, then g'=ae w and f =2a cos 0. That is, if P 
and if" are corresponding points in the f and z' planes, 
respectively, then as P traverses the z axis from 2a to 
—2a, P' traverses the circle ae i$ from 0=0 to 0=ir } 
and as P retraces its path to f =2a, P' completes the 
circle. The transformation (5) then may be seen to 
map the entire £ plane external to the line 4a uniquely 
into the region external (or internal) to the circle of 
radius a about the origin in the z' plane. 

Let us invert equations (6) and solve for the elliptic 
coordinates if and 0. (Fig. 3.) We have 




/ Plane z' Plane 

Floras 3.— Transformation by elliptic coordinates 



cosh ^= 



2a cos 0 

sinh 

and since cosh V>-sinh 2 ^= 1 

( * Y-7 y Y-i 

\2a cos 0/ \2a sin 0/ 
or solving for sin 2 0 (which can not become negative), 



2a sin 0 



where 
Similarly we obtain 



2sin 2 0=y+ A / P 2 +(|y 

' --(£)'-(£)' 



(7) 



( * \\f y V«i 

\2a cosh ^/ \2a smh 



v,2a cosh ^. 
or solving for sinh V 
2 sinh V 



(8) 



We note that the system of radial lines 0 = constant 
become confocal hyperbolas in the f plane. The circles 
f = constant become ellipses in the f plane with major 
axis 2a cosh # and minor axis 2a sinh f . These orthog- 
onal systems of curves represent the potential lines and 
streamlines in the two planes. The foci of these two 
confocal systems are located at (± 2a, 0). 

408318 0—41 2 



Equation (8) yields two values of.^ for a given 
point (x, y), and one set of these values refers to the 
correspondence of (x, y) to the point (ae*, 0) external to 
a curve and the other set to the correspondence of 
(x, y) to the point (ae - *, — 0) internal to another curve. 
Thus, in Figure 3, for every point external to the 
ellipse Ei there is a corresponding point external to the 
circle C\, and also one internal to O x '. 

The radius of curvature of the ellipse at the end of 

the major axis is p = 2a ^ or for small values of if/, 

ps^a^ 2 . The leading edge is at 

2a cosh ^sg2a^l + y )s*2a. + ~ 

Now if there is given an airfoil in the £ plane (fig. 4), 
and it is desired to transform the airfoil profile into a 
curve as nearly circular as possible in the z' plane by 
using only transformation (5), it is clear that the axes 
of coordinates should be chosen so that the airfoil 
appears as nearly elliptical as possible with respect to 
the chosen axes. It was seen that a focus of an 
elongated ellipse very nearly bisects the line joining 
the end of the major axis and the center of curvature 
of this point; thus, we arrive at a convenient choice of 
origin for the airfoil as the point bisecting the line of 
length 4a, which extends from the point midway be- 
tween the leading edge and the center of curvature of 
the leading edge to a point midway between the 
center of curvature of the trailing edge and the trailing 
edge. This latter point practically coincides with the 
trailing edge. 

The curve B, defined by ae* +u , resulting in the z ! 
plane, and the inverse and reflected curve B', defined 
by ae~*~ w , are shown superposed on the f plane in 
Figure 4. The convenience and usefulness of trans- 




Figure 4. — Transformation of airfoil into a nearly circular contour 



formation (5) and the choice of axes of coordinates 
will become evident after our next transformation. 



The transformation s' = ze 



.—Consider the trans- 



formation z' = ze /ll) where/(z) — 2 Each exponential 

o * 

Cn 

term e ^represents the uniformly convergent series 



i 
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where the coefficients c n =A„ + iB n are complex num- 
bers. For /(s) convergent at all points in a region 
external to a certain circle, z' has a unique real abso- 
lute value [z|e |/M l in the region and its imaginary part 
is definitely defined except for integral multiples of 
2ir£. When 2= <» , z' = ze e °. The constant Co — A a + 
iB 0 is then the determining factor at infinity, for the 
field at infinity is magnified by e 4 ° and rotated by the 
angle B 0 . It is thus clear that if it is desired that the 
regions at infinity be identical, that is, 2' = z at infinity, 
the constant c 0 must be zero. The constants c t and e t 
also play important roles, as will be shown later. 

We shall now transform the closed curved 7 z' = ae* +i ° 
into the circle 2=ae* 0+1 * (radius ae*°, origin at center) 
by means of the general transformation (reference 2) 



1 2" 



(10) 



■ze 



which leaves the fields at infinity unaltered, and we 
shall obtain expressions for the constants A n , B„, and 
i^o. The justification of the solution will be assured by 

the actual convergence of 2^?, since if the solution 

1 2" 

exists it is unique. 

By definition, for the correspondence of the bound- 
ary points, we have 

a *-*o+i(e- v ) (10') 



z —ze 



Also 

Consequently 



. S(A+iB.)^ 
1 • 



where 



z=ae* 0 +tv 



On writing 2=B(cos <p+i sin <p) where R — ae* a , we 
have 

" 1 

4>— ^0 + i(fi — <p) = 2 (An + iB„)-^(coa n<p — i sin n<p) 

Equating the real and imaginary parts of this relation, 
we obtain the two conjugate Fourier expansions: 

tf-*o=S cos nv+j* sin n?] (11) 



■ VB A 1 

8-<p=2\j£ cos «*>-g5 8"^ ™*»J 



(12) 



From equation (11), the values of the coefficients jfi> 
jpt> and the constant f 0 are obtained as follows: 



"ik=- cos n<P dtp 
u * 0 



5i=- fty sin n<p dtp 



(a) 
0>) 



7 Unless otherwise stated, and 9 will now be used in this restricted sense, i, e., as 
defining the boundary curve itself, and not all points in the i' plane, 



1 2ir 
^o=2^./Vd<9 



(c) 



The evaluation of the infinite number of constants 
as represented by equations (a) and (b) can be made 
to depend upon an important single equation, which 
we shall obtain by eliminating these constants from 
equation (12). 

Substitution of (a) and (b) for the coefficients of 
equation (i2) gives 

2» [~ 2x 

(0— <p)' =- 2 cos n<p' f $(<p) sin n<p dtp 

r 1 L o 

2, 1 

— sin n<p' f $(<p) cos n<p dtp 

0 J 

where ^(<p)—^ and (fi—tp)' represents 0 — tp as a func- 
tion of <p' , and where tp' is used to distinguish the angle 
kept constant while the integrations are performed. 
The expression may be readily rewritten as 



($—<p)'=- 2 ,/>(?>) (sin n<p coa n<p f — co9ntpsian<p')dtp 
i <■> 2ir 



But 



=- 2 f f(tp) sin n(<p— <p') dtp 
* 1 0 



<p-<P 



COS (271+ 1) 



2 sin ra (*>-*>') =2 cot 2 



sm 



Then 



(*-*)'■ 



lim 
, 2 tt 



2t 



£ S*to) cot Z-^- d v 

cos(2n+l) ^fi 

- — d v 



2-k 



SHv) 

0 



sm 



<P~<P 



The first integral is independent of n, while the latter 
one becomes identically zero. 

Then finally, representing <p-d by a single quantity 
£, viz <p—d=e=e(<p), we have 

2r 



e(«>') = -AStiv) cot 



(13) 



By solving for the coefficients in equation (12) and 
substituting these in equation (11) it may be seen that 
a similar relation to equation (13) holds, for the func- 
tion i(<p). 



2t 



2t 



The last term is merely the constant fa, which is, as 
has been shown, determined by the condition of mag- 
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nification of the s and s' fields at infinity. The 



I 



2x 



corresponding integral s - f t{<p) d<p does not appear in 



0 



equation (13), being zero as a necessary consequence 
of the coincidence of directions at infinity and, in 
general, if the region at infinity is rotated, is a constant 
different from zero. 

Investigation of equation (13). — This equation is 
of fundamental importance. A discussion of some of 
its properties is therefore of interest. It should be 
first noted that when the function i>(<p) is considered 
known, the equation reduces to a definite integral. 
The function 8 e(p) obtained by this evaluation is the 
"conjugate" function to so called because of the 
relations existing between the coefficients of the 
Fourier expansions as given by equations (11) and (12). 
For the existence of the integral it is only necessary 
that be piecewise continuous and differentiable, 
and may even have infinities which must be below 
first order. We shall, however, be interested only in 
continuous single-valued functions having a period 2t, 
of a type which result from continuous closed curves 
with a proper choice of origin. 

If equation (13) is regarded as a definite integral, it is seen 
to be related to the well-known Poisson integral which solves 
the following boundary- value problem of the circle. (Reference 
3.) Given, say for the z plane a single-valued function u(B,t) 
for points on the circumference of a circle w=Re tr (center at 
origin), then the single-valued continuous potential function 
u(r,a) in the external region 2=re ir of the circle which assumes 
the values u (R, t) on the circumference is given by 

1 i^—R* 
uM =2^ u(Rir) rW-2Kr cos (<r-r) d * 

and similarly for the conjugate function v(r,<r) 



i 2t 



T*-R? 



0 



B 2 +r'—2Br cos (o--t) 



dr 



These may be written as a single equation 

u(r,<r) +iv(r,<r) -/(«) =i £/(«) S±H«J« 

where the value f{z) at a point of the external region «=^re ,v is 
expressed in terms of the known values /(to) along the circum- 
ference v>=Re iT . In particular, we may note that at the 



boundary itself, since »' g<r fr ■ 
2* 



cot 



^ i we have 



«(fl,<r) +iv(R,<r) - MR,T) +iv(R,r)] Cot ^—^ dr > 



which is a special form of equations (13) and (14). 

The quantity ^ is immediately given as a function 
of 0 when a particular closed curve is preassigned, and 
this is our starting point in the direct process of trans- 
forming from airfoil to circle. We desire, then, to find 
the quantity ^ as a function of <f> from equation (13), 
and this equation is no longer a definite integral but an 

« This function will be called "confnrmal angular distortion" function, for reasons 
evident later. 



integral equation whose process of solution becomes 
more intricate. It would b& surprising, indeed, if 
anything less than a functional or integral equation 
were involved in the solution of the general problem 
stated. The evaluation of the solution of equation (13) 
is readily accomplished by a powerful method of suc- 
cessive approximations. It will be seen that the 
nearness of the curve ae**** to a circle is very signifi- 
cant, and in practice, for airfoil shapes, one or at most 
two steps in the process is found to be sufficient for 
great accuracy. 

The quantities 4> and e considered as functions of <p 
have been denoted by and e(<p), respectively. 
When these quantities are thought of as functions of 8 
they shall be written as vK0) and e(fi), respectively. 

Then, by definition 



and 



He) ««[*(*)] 

Since 0=e, we have 

<p(8) = 0-+i(0) J 



(15) 



(16) 



We are seeking then two functions, $(<p) and «(<?), 
conjugate in the sense that their Fourier series expan- 
sions are given by (11) and (12), such that \fr[(»(0)] = 
i(8) where ?(0) is a knovm single-valued function of 
period 2t. 

Integrating equation (13) by parts, we have 



0 

<P—<p' 



The term log sin g is real only in the range <p = <p' to 

¥> = 27r+p', but we may use the interval 0 to 2x for <p 
with the understanding that only the real part of the 
logarithm is retained. 
Let us write down the following identity: 

. . <p—ip' , . 8-9' 
log sin 2L 2 - JC - =log sin — g— 



+log- 



sin 



(e+i I )-(9+« 1 )' 



sin 



8-8' 



-+log 



sin- lJ —z± iJ - 

(8+- ei )-(d+W 



sin 



sin 



(9 + « t )-(0 + e>)' 



• +log 



■ (9 + e>-l)-(9 + i*-i) 



> + 



sin 



(0-^-e^,)-(0 + i B ) , 



sin 



(17) 



(0-f e)-(0 + e)' 



+ log 



(0+i a . i )-(0+^,. 1 ) 
sin 2 



7+log 



where in the last term we recall that 0 + «(0) = <p(9) ; and 
where it may be noted that each denominator is the 
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numerator of the preceding term. The symbols e k 
(k = l, 2, • • •, n) represent functions of 0, which thus 
far are arbitrary. 9 

Since by equation (15) 7(0) 1B tf?(9)] w e have for 
corresponding elements d0 and dtp 

Then multiplying the left side of equation (17) by 
~ d*> and the right side by ^ d0 and inte- 

grating over the period 0 to 2ir we obtain 



i[#)]-W = i/log3in^ d -f d»+ . . 



2t 



sin 



d^W 



2 - 3in (0+«(0))-(0+i(0))' 

+ -t 0 ;log ^5)3±r d -^ de (18) 



where 4 = 1, 2, .. . ., n. 

We now choose the arbitrary functions l k (0') so that 

i„(0')=O 

and 

«>(»') =\ S log ^ »+ VJ- g+^J' jffi dfl (19) 

where 2, . . ., ra. 
Equation (18) may then be written 

i(0')=«b + *i+(«2-*i) • • • +6.-*,-i) + (?-*») (20T 
or l(0') = Xi + Xs + . . 

where X*(0') = e*_i and is in fact the k th term of 
equation (18). The last term we denote by X. 

From equation (19) we see that the function «*(#') is 
obtained by a knowledge of the preceding function 
e*_i(0')- For convenience in the evaluation of these 
functions, say 

wr>-i/io. sm^tM^r dg) d , 



d0 



we introduce a new variable <p k denned by 

?*(0)=0+**(0) (t-i, 2, . . ., a) 

Then 



-i /log sin 

From the definition of <p„ as 
*>*(«) 



d<p* 



dip* 



(21) 



• The symbol (»+*»)' represents and is used to denote the same function 

of ef that «+•*(« is of 9. The variables 0 and C are regarded as independent of each 
other. 



we may also define the symbol e k (<p k ) by 
where 

-«*[**(*)] 

It is important to note that the symbols t k , e*, «** 
denote the same quantity considered, however, as a 
function of 0, <p k , <p k - U respectively. 

The quantities («» — e*_i) in equation (20) rapidly 
approach zero for wide classes of initial curves i£(0), 
i. e., ${B(tp k )] very nearly equals ?[0(p*-h)] for even 
small k's. The process of solution of our problem is 
then one of obtaining successively the functions i£(0), 
ROM], 'ffiffa)], • • • • where $[0(<p n )] and 

*J0(<Pn)] become more and more "conjugate." The 
process of obtaining the successive conjugates in prac- 
tice is explained in a later paragraph. We first pause 
to state the conditions which the functions <p k are sub- 
ject to, necessary for a one-to-one correspondence of 
the boundary points, and for a one-to-one corre- 
spondence of points of the external regions, i. e., the 
conditions which are necessary in order that the 
transformations be conformal. 

In order that the correspondence between boundary 
points of the circle in the z plane and boundary points 
of the contour in the z' plane be one-to-one, it is 
necessary that 8(<p) be a mono tonic increasing function 
of its argument. This statement requires a word of 
explanation. We consider only values of the angles 
between 0 and 2ir. For a point of the circle boundary, 
that is, for one value of <p there can be only one value 
of 9, i. e., 6{<p) is always single valued. However, ^(0), 
in general, does not need to be, as for example, by a 
poor choice of origin it may be many valued, a radius 
vector from the origin intersecting the boundary more 
than once; but since we have already postulated that 
f(0) is single valued this case can not occur, and <p(6) 
is also single valued. If we decide on a definite direc- 
tion of rotation, then the inequality g~ ^ 0 expresses 

the statement that as the radius vector from the origin 
sweeps over the boundary of the circle C, the radius 
vector in the z' plane sweeps over the boundary of B 
and never retraces its path. 
The inequality 



corresponds to 



*L =1 _4i(£)^ 

dtp dtp 



Also, the condition 



dp 



<1 



d0 1 + 



de - 



corresponds to 



de(0) 
d0 - 
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Multiplying by we get 



This relation is shown in Figure 5 as a rectangular 




Figure 5.— The quantity as a function of 

hyperbola. We may notice then that the monotonic 
behavior of <p{6) and 8(<p) requires that remain on 
the lawer branch 10 of the hyperbola, i. e., 

(22) 

It will be seen later that the limiting values 



dip 



, de(yj) /. di de A 



correspond to points of infinite velocity and of zero 
velocity, respectively, arising from sharp corners in the 
original curve. 

The condition for a one-to-one conformal corre- 
spondence between points of the external region of the 
circle and of the external region of the contour in the 
z' plane may be given (reference 5, p. 98 and reference 
6, Part II) as follows: There must be a one-to-one 
boundary point correspondence and the derivative of 

Z™ 

the analytic function z' = ze l given by equation (10) 
must not vanish in the region. That is, writing g(z) 

for 2^ we have 

1 Z* 

g^=e' (i) (l + z^) 7*0 for \z\>R or since 

the integral transcendental function e"- l) does not vanish 
in the entire plane, the condition is equivalent to 
_ dg(g) 



dz 



'-9*-l for \z\>R 



'« The values of the upper branch of the hyperbola arise when the region internal 
to the curve at *«• is transformed into the external region of a circle, but may also 
there be avoided by defining <=»>+» instead of «>-0. 



By equation (10') we have on the boundary of the 
circle, g(Re iv ) = - fa - it., and 

dg(z) M d[»(ri-M(»)] 
2 dz ~ m iRe'* d 9 

_ dt (f) j d t(f>) 
d<p dip 

the first term on the right-hand side being real and the 
last term a pure imaginary. We have already postu- 
lated the condition 

-de 

— 00 a t~ = i 
dip 

as necessary for a one-to-one boundary point corre- 
spondence. Now by writing z = £ + mj and g ^^ °= 
P (£>*/) +iQ(Z,v), we note that gives the boundary 

values of a harmonic function P (£,»?) and therefore this 
function assumes its maximum and minimum values 
on the boundary of the circle itself. (Reference 3, p. 



Hence 



223.) 
external region, i 



dg(z) 

z dT can 



e., 



dz' 
dz 



never become — 1 in the 
can never vanish in this 



region. 

At each step in the process of obtaining the succes- 
sive conjugates we desire to maintain a one-to-one 
correspondence between 8 and ip k , i. e., the functions 
0(*>*) and ipic(0) should be monotonic increasing and are 
hence subject to a restriction similar to equation (22), 
viz, 

„-=rd«* 



'dip, 



<1 



(22') 



The process may be summed up as follows : We con- 
sider the function ij/(6) as known, where <K0) is the 
functional relation between ^ and 0 defining a closed 
curve ae* +tf . The conjugate of $(8) with respect to 8 
is ii(0). We form the variable ipi = 8 + ti(6) and also 
the function i£[0(<pi)]. The conjugate of i£[0(«5i)] with 
respect to ipi is e*2(<Pi) which expressed as a function of 
8 is i 2 (8). We form the variable (o 2 = 0 + 6 2 (0) and the 
function ^[8(ip 2 )]. The conjugate of $[8(ip 2 )] is. t* 3 (ip 2 ), 
which as a function of 8 is e 3 (0), etc. The graphical 
criterion for convergence is, of course, reached when 
the function f[0 (ip„)] is no longer altered by the 
process. The following figures illustrate the method 
and exhibit vividly the rapidity of convergence. The 
numerical calculations of the various conjugates are 
obtained from formula I of the appendix. 

In Figure 6, the $(0) curve represents a circle re- 
ferred to an origin which bisects a radius (obtained 
from an extremely thick Joukowsky airfoil) (see p. 26) 
and has numerical values approximately five times 
greater than occur for common airfoils. The ^/(ip) 
curve is known independently and is represented by 
the dashed curve. The process of going from $(0) to 
4>{ip) assuming if>(ip) as unknown is as follows: The 
function «i(0), the conjugate function of $(0), is found. 
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The quantity ^ is then plotted against the new variable 
<Pi = 8 + e 1 (d) (i. e., each point of \j/(d) is displaced hori- 
zontally a distance e0 and yields the curve jfr [0(*>i)]. 
(likewise, ^(0) is plotted against <pi yielding «i ($*).) 



is drawn at P'. This process yields the function ^ (ft) . 
The quantity $ is now plotted against the new variable 
p 2 = 0+e 2 (0) (i. e., each point of ${d) is displaced hori- 
zontally a distance e 2 ) giving the function $[0(e»)]. 




2.5 3.0 W-jt 3.S 4.0 
Argument (0, <p,,<p,,<p in radians) 



6.0 2k 



The function «* 2 (#>i) is now determined as the conjugate 
function of #[0(*>i)]. This function expressed as a 
function of 0 is «* 2 [v»i (0)] =^(0) . It is plotted as follows : 



Fionas 6.— The process of obtaining successive conjugates 

This curve is shown with small circles and coincides 
with $(<p). Further application of the process can 
yield no change in this curve. It may be remarked 
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Fiqvbk 7.— Process applied to transforming a square into a circle 



2n 



At a point P of e* 2 (^>i) and Q of fifa) corresponding to 
a definite value of <pi one finds the value of 0 which 
corresponds to <p t by a horizontal line through Q meet- 
ing ti(fi) in Q'; for this value of 0, the quantity e 2 at P 



here- that for nearly all airfoils used in practice one 
step in the process is sufficient for very accurate results. 

As another example we shall show how a square 
(origin at center) is transformed into a circle by the 
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method. In Figure 7 the #(0) curve is shown, and in 
Figure 8 it is reproduced for one octant. 11 The value 
is #(0)=log sec 0. The function f[d(<f>i)] is shown 
dashed; the function $[0(tt>)] is shown with small 
crosses; and ^[flfe)] is shown with small circles. The 
solution ^(<p) is represented by the curve with small 
triangles and is obtained independently by the known 
transformation (reference 3, p. 375) which transforms 
the external region of a square into the external region 
of the unit circle, as follows: 



W (g) = /VEEIr 



26 



2 s 



where jpQ) denotes a power series. Comparing this 

with equation (10), we find that 4>{<p) except for the 

constant & is given as the real part of log^l+P^-^J 

evaluated for z=e i *, and that «(«>) is given as the 
negative of the imaginary part. It may be observed 
in Figure 8 that the function £[0(?a)] very nearly 
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Fiqube 8.— Process applied to transforming a square into a circle 

equals The functions «(<e) and i(0) are shown in 

Figure 7 (a); we may note that at ¥>=|, which corre- 
sponds to a comer of the square, ^=1 or also, 
lie * 

ar-- 

" Because of the symmetry involved only the interval 0 to \ need be used. The 
integral in the appendix can be treated as 



«(*>0 — £S*M cot 



2*{"™ 2 



<1<P 



"~*0 Wril^'^-^-^^HVJld* 



It may be remarked that the rapidity of convergence 
is influenced by certain factors. It is noticeably af- 
fected by the initial choice of eb(0). The choice 
eo(0) =0 implies that 0 and <p are considered to be very 
nearly equal, i. e., that oe*" 1 " 49 represents a nearly cir- 
cular curve. The initial transformation given by 
equation (5) and the choice of axes and origin were 
adapted for the purpose of obtaining a nearly circular 

S, 




P- 




(b) 



Figure 9.— Translation by the distance OM 

curve for airfoil shapes. If we should be concerned 
with other classes of contours, more appropriate 
initial transformations can be developed. If, how- 
ever, for a curve ae*" 1 "* 9 the quantity t = <p— 0 has large 
values, either because of a poor initial transformation 
or because of an unfavorable choice of origin, it may 
occur that the choice «o(0) = O will yield a function 

ei(<Pi) for which -j^ 5 - may exceed unity at some points, 
a<pi 

thus violating condition (22'). Such slopes can be 
replaced by slopes less than unity, the resulting func- 
tion chosen as 4^(0) and the process continued as 
before. 12 Indeed, the closer the choice of the function 
e\)(0) is to the final solution i(0), the more rapid is the 
convergence. The case of the square illustrates that 
even the relatively poor choice «o(0) = 0 does not appre- 
ciably defer the convergence. 

The translation z t = z+ci. — Let us divert our 
attention momentarily to another transformation 
which will prove useful. We recaE that the initial 
transformation (eq. (5)) applied to an airfoil in the f 
plane gives a curve B in the z' plane shown schemati- 
cally in Figure 9(a). Equation (10) transforms this 
curve into a circle C about the origin 0 as center and 
yields in fact small values of the quantity <p—9. We 
are, however, in a position to introduce a convenient 
transformation, namely, to translate the circle O into 
a most favorable position with respect to the curve B 
(or vice versa). These qualitative remarks admit of a 
mathematical formulation. It is clear that if the 
curve B itself happens to be a circle 13 the vector by 
which the circle O should be translated is exactly the 
distance between centers. It is readily shown that 

u The first step in the process is now to define va-e+mW) and form the function 
?[*(**>)]• The conjugate function of is «•«{*») which expressed as a function 

of I is «(«). etc. . 

» See p. 28. 



14 



BEPOBT NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



then equation (10) should contain no constant term. 
We have 



2 =26 



1 z* 



(10) 



(l+^+ • .)etc. 
= Z (l+|+|+ ) {10a) 



where 14 



&i = d 

*2 = C 2 + y 

7 , I C l 3 

#3 = C 3 + CgCi + "g 



It is thus apparent that if equation (10) contains no 
first harmonic term, i. e., if 

Ci^Ai+iBi — -f4ie'*d<p = 0, 
1 0 

the transformation is obtained in the so-called normal 
form 



2, Zi' 



(23) 



This translation can be effected either by substituting 
a new variable Z\ — z+ c u or a new variable z x ' = 2' — c x . 



10 
.05 
£ 0 
-.05 
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Figure 10.— The £(«) and *,(«,) carves (for Clark Y airfoil) 

This latter substitution will be more convenient at 
this time. Writing 



2,'=ae*»+«fc c l = aei'+", and z' = ae*+ if 



we have 



ae 4,i+tm = ae ++u_ aey +ui 



The variables and (7,, can be expressed in terms of 
0, y, and 5. Tn Figure 9(b), P is a point on the B 

» These constants can be obtained in a recursion form. See footnote 18, 



curve, i. e., OP = ae*, PQ represents the translation 
vector Ci = aet+ is , OQ is ae* l+i>l , and angle POQ is 
denoted by Then by the law of cosines 



e ».-e 2 * + « 2 T-2eV cos (0-5) 
and by the law of sines 

= e sin (0-5) 
sin jt ^ 



(a) 



or 



0! = 0+ju = 0+tan 



_ t sin (0-8) 



(b) 



C os (0-5) 

In Figure 10 are shown the ?(0)^and e(0) curves for the 
Clark Y airfoil (shown in. fig. 4) and the fc(6\) and 
«i(0i) curves which result when the origin is moved 
from 0 to M. It may be noted that e^) is indeed 
considerably smaller than e(0). It is obtained from 



(*-W--s- f M<p) cot 

and the constant ^ 0 is given 15 by 
1 2* 



<P~<P 



0 



The combined transformations. — It will be useful to 
combine the various transformations into one. We 
obtain from equations (5) and (10) an expression as 
follows: 



r=2d cosh ^log 2+ 2 



(24) 



or we can also obtain a power series development in z 

? -r - • ' (25) 

where 19 a n = k n+i + a 2 h„-i 



2 

On^^+l" 

The constants k„ may be obtained in a convenient 



recursion form as 



*l = Ci 

2ki = kiCi + 2c2 

3& 3 = A 2 Ci + 2^02 + 3c 3 

4A 4 = kiCi + 2k^ 2 + 3kiC 3 + 4c 4 



The constants h n have the same form as k n but with 
each C( replaced by — c t (and Ao = l). It will be re- 

•» The constant *h is invariant to change of origin. (See p. 26.) It should be 
remarked that the translation by the vector a is only a matter of convenience and 
is especially useful for very irregular shapes. For a study of the properties of airfoil 
shapes we shall use only the original t(v) curve. (Fig. 10(a).) 
By'equations (5) and (10) we have 

m t. 00 n 

s — —2 — 

• 

1 i 2 n 

The constant*, is thus the coefficient of — in the expansion of e 1 andthe constant 



An the coefficient of — in the expansion of e 



l 2 



For the recursion form for *„ 
see Smithsonian Mathematical Formulas and Tables of Elliptic Functions, p. 120. 
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called that the values of c„ are given by the coefficients 
of the Fourier expansion of \f>(<p) as 



and 



c 1 2t 

ij? = i / ^(^e^Av where R = ae*« 



The first few terms of equation (25) are then as 
follows : 

ci 2 «,» 



, c 2 + 2 +a 2 , c 3 + c 2 Ci+ 6 

f-2 + C, + + . J 



(25') 



By writing s l = z + ci, equation (25) is cast into the 
normal form 



(26) 



The constants 6„ may be evaluated directly in terms 
o.f a„ or may be obtained merely by replacing tiv) by 
(<p) in the foregoing values for a n . 

The series given by equations (25) and (26) may be 
inverted and z or z x developed as a power series in f . 
Then 

2(0 =f-Ci--p p p . . . (27) 

and 



ZlUJ - S f _ f 2 i" 3 



(28) 



The various transformations have been performed 
for the purpose of transforming the flow pattern of a 



\ 



I 



.... \ -zii^^ 


>^ 









Figure U.— Streamlines about circle with zero circulation (shown by the full 
lines) Q= - Vsinh m sin v= constant 

circle into the flow pattern of an airfoil. We are thus 
led immediately to the well-known problem of deter- 
mining the most general type of irrotational flow 
around a circle satisfying certain specified boundary 
conditions. 

The flow about a circle. — The boundary conditions 
to be satisfied are: The circle must be a streamline of 
flow and, at infinity, the velocity must have a given 
magnitude and direction. Let us choose the £ axis as 
corresponding to the direction of the velocity at 

408318 O— 41 3' 



infinity. Then the problem stated is equivalent to 
that of an infinite circular cylinder moving parallel to 
the | axis with velocity V'm a fluid at rest at infinity. 

The general complex flow potential 17 for a circle of 
radius B, and velocity at infinity V parallel to the x 
axis is 

w(2) = -v( s+ |)-glog| _ (29) 

where T is a real constant parameter, known as the 




Figure 12.— Streamlines about circle for V=o Q=— gy»=constant 

circulation. It is defined as cjw,ds along any closed 
curve inclosing the cylinder, v, being the velocity 
along the tangent at each point. 

"Writing z = Re» +i * and w=P+iQ, equation (29) be- 
comes 



or 



w = — V cosh(/i + itp) — 2^(m + i(p) 
P= — V cosh fx cos <P + ^<P 

r 

Q= — V sinh n sin •f—^v- 

For the velocity components, we have 
dw . T/ /, R 2 \ 



(29') 



ir 

2wz 



(30) 



In Figures 11 and 12 are shown the streamlines for 
the cases r = 0, and V=Q, respectively. The cylinder 
experiences no resultant force in these cases since all 
streamlines are symmetrical with respect to it. 

The stagnation points, that is, points for which u 

and v are both zero, are obtained as the roots of -gj = 0. 
This equation has two roots. 

~iiv 

and we may distinguish different types of flow accord- 
ing as the discriminant le^i? 2 ^ 2 — r 2 is positive, zero, 
or negative. We recall here that a conformal trans- 
formation w =/(s) ceases to be conformal at points 

where vanishes, and at a stagnation point the flow 

divides and the streamline possesses a singularity. 

" Reference 4, p. 56 or reference 5, p. 118. The log term must be added because 
the region outside the infinite cylinder (the point at infinity excluded) is doubly 
connected and therefore we must include the possibility of cyclic motion. 
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The different types of flow that result according 

as the parameter r 2 ^ lQnZRW 2 are represented in 

Figure 13. In the first case (fig. 13 (a)), which will not 
interest us later, the stagnation point occurs as a 
double point in the fluid on the ij axis, and all fluid 
within this streamline circulates in closed orbits around 
the circle, while the rest of the fluid passes downstream. 
In the second case (fig. 13 (b)), the stagnation points 
are together at S on the circle Re i<p and in the third 
case (fig. 13 (c)) they are symmetrically located on the 
circle. We have noted then that as r increases from 
0 to 4tBV the stagnation points move downward on 

the circle Re t<p from the £ 



axis toward the ij axis. 
Upon further increase in 
r they leave the circle and 
are located on the n axis in 
the fluid. 

Conversely, it is clear 
that the position of the 
stagnation points can de- 
termine the circulation r. 
This fact will be shown to 
be significant for wing- 
section theory. At pres- 
ent, we note that when 
both r and 0 a marked 
dissymmetry exists in the 
streamlines with respect to 
the circle. They are sym- 
metrical about the ij axis 
but are not symmetrical 
about the £ axis. Since 
they are closer together on 
the upper side of the circle 
than on the lower side, a 
FiouBE is.-streamiines about circle resultant force exists per- 

[from Lagally— Handbuoh der Physik pendicular to the motion. 

Bd.virj Q=T'sinh/isin t .-^i-con. We shall now combine 
stant(a)r«>i6i»B'V>(b)rj-i6iijpvi the transformation (27) 
M n< iftwv. ^ the flow formula for 

the circle equation (29) and obtain the general complex 
flow potential giving the 2-dimensional irrotational flow 
about an airfoil shape, and indeed, about any closed 
curve for which the Riemann theorem applies. 

The flow around the airfoil. — In Figure 14 are 
given, in a convenient way, the different complex 
planes and transformations used thus far. The com- 
plex flow potential in the z plane for a circle of radius 
R origin at the center has been given as 




W ( 2 )--y( 3+ D-|£io gi 



(29) 



permit of a change in the direction of flow at infinity 
by the angle a which will be designated angle of attack 
and defined by the direction of flow at infinity with 
respect to a fixed axis on the body, in this case the 
axis v>=0. This flow is obtained simply by writing 
ze ta for z in equation (29) and represents a rotation of 



S Plane 



where V, the velocity at infinity, is in the direction of 
the negative £ axis. Let us introduce a parameter to 




Floras 14.— The collected transformations 

the entire flow field about the~circle by angle a. We 
have 

w(z) = -v(ze i "+ E ^e-*'^- i £logz (31) 
dw 



dz 



Ve V z* e J 2, 



%-KZ 



(32) 



Since a conformal transformation maps streamlines 
and potential lines into streamlines and potential lines, 
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we may obtain the complex flow potentials in the 
various planes by substitutions. For the flow about 
the circle in the z x plane, z is replaced by z x — c x 



w{zi) 



-A 



(2, -<?,)**■+ 



]-glog(*- 



Cl ) (31') 



(21 -ci). 

dF, TV-fl ^"" l ^— (32') 

For the flow about the B curve in the 2' plane, z is 
.replaced by z(z') (the inverse of eq. (10a)) and for the 
flow about the airfoil in the f plane z is replaced by 
2(f) from equation (27) 

W® = - V[«(r)«<-+ J|ye-'«] log 2 (f) (33) 



242(f)] 



(34) 



The flow fields at infinity for all these transformations 
have been made to coincide in magnitude and direction. 

At this point attention is directed to two important 
facts. First, in the previous analysis the original 
closed curve may differ from an airfoil shape. The 
formulas, when convergent, are applicable to any 
closed curve satisfying the general requirements of 
the Riemann theorem. However, the peculiar ease of 
numerical evaluations for streamline shapes is note- 
worthy and significant. The second important fact is 
that the parameter T which as yet is completely unde- 
termined is readily determined for airfoils and to a 
discussion of this statement the next section is devoted. 
It will be seen that airfoils may be regarded as fixing 
their own circulation. 

Kutta-Joukowsky method for fixing the circula- 
tion. — All contours used in practice as airfoil profiles 
possess the common property of terminating in either 
a cusp or sharp corner at the trailing edge (a point of 
two tangents). Upon transforming the circle into an 

Idz 

airfoil by f =f(z), we shall find that ^ is infinite at 

the trailing edge if the tail is perfectly sharp (or very 
large if the tail is almost sharp). This implies that 



the numerical value of the velocity 



dw 
dz 



--\v\ is 



dw 
dz 



infinite (or extremely large) provided the factor 

is not zero at the tail. There is but one value of the 
circulation that avoids infinite velocities or gradients 
of pressure at the tail and this fact gives a practical 
basis for fixing the circulation. 

The concept of the ideal fluid in irrotational poten- 
tial flow implies no dissipation of energy, however large 
the velocity at any point. The circulation being a 
measure of the energy in a fluid is unaltered and inde- 
pendent of time. In particular, if the circulation is 
zero to begin with, it can never be different from zero. 



infinite. It is clear that the factor 



must then be 



However, since all real fluids have viscosity, a better 
physical concept of the ideal fluid is to endow the 
fluid with infinitesimal viscosity so that there is then 
no dissipation of energy for finite velocities and pres- 
sure gradients, but for infinite velocities, energy losses 
would result. Moreover, by Bernoulli's principle the 
pressure would become infinitely negative, whereas a 
real fluid can not sustain absolute negative pressures 
and the assumption of ^compressibility becomes in- 
valid long before this condition is reached. It should 
then be postulated that nowhere in the ideal fluid from 
the physical concept should the velocity become 

dw 
dz 

zero at the trailing edge in order to avoid infinite 
velocities. It is then precisely the sharpness of the 
trailing edge which furnishes us the following basis for 
fixing the circulation. 

It will be recalled that the equation gj = 0 deter- 
mines two stagnation points symmetrically located on 
the circle, the position of which varies with the value 
of the circulation and conversely the position of a 
stagnation point determines the circulation. In this 
paper the x axis of the airfoil has been chosen so that 
the negative end (0 = jt) passes through the trailing 
edge. From the calculation of *=><p—0 (by eq. (13)) 
the value of ip corresponding to any value of 8 is deter- 
mined as ip = 6 + «, in particular at 0 = ?r, <p — ir + 0, where 
/3 is the value of « at the tail and for a given airfoil is a 
geometric constant (although numerically it varies 
with the choice of axes). This angle /S is of consid- 
erable significance and for good reasons is called the 
angle of zero lift. The substance of the foregoing 
discussion indicates that the point z*=Be t(T+ & = —Re i0 
is a stagnation point on the circle. Then for this value 
of 2, we have by equation (32) 



2irz 



= 0 



or 



r= -27^i^F^'e ^( " + < i, (1 - e'^+f") 



=4t#v(- 



2i 



= 4irflFsin (a + j3) 



(35) 



This value of the circulation is then sufficient to 
make the trailing edge a stagnation point for any value 
of a. The airfoil may be considered to equip itself 
with that amount of circulation which enables the 
fluid to flow past the airfoil with a minimum energy 
loss, just as electricity flowing in a flat plate will dis- 
tribute itself so that the heat loss is a minimum. The 
final justification for the Kutta assumption is not only 
its plausibility, but also the comparatively good agree- 
ment with experimental results. Figure 15 (b) shows 
the streamlines around an airfoil for a flow satisfying 
the Kutta condition, and Figures 15 (a) and 15 (c) illus- 
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trate cases for which the circulation is respectively too 
small and too large, the stagnation point being then on 
the upper and lower surfaces, re- 
spectively. For these latter cases, 
the complete flow is determinable 
only if, together with the angle of 
attack, the circulation or a stag- 
: '(b) nation point is specified. 

Velocity at the surface. — The 
flow formulas for the entire field 
are now uniquely determined by 
substituting the value of V in equar 

noun i 6 .-(a) now with tions (33) and . (34) ■ We are, how- 
circuiation smaller than for ever, in a position to obtain much 

Kutta condition; (b) flow „•_>„_ . j „ ■ , 

satistyingKuttaoondition; simpler and more convement re- 
(c) flow with circulation lations for the boundary curves 
SSSSm*"" fM KU " a themselves. Indeed, we are chiefly 
interested in the velocity at the 
surface of the airfoil, which velocity is tangential to 
the surface, since the airfoil contour is a streamline of 
flow. The numerical value of the velocity at the 
surface of the airfoil is 




did 


dw 




dz 




dz' 




dz 




dz' 




df 



We shall evaluate each of these factors in turn. From 
equations (32 and (35) 

dz Ve V s 2 ) 2wz 
At the boundary surface z=Re i,f , and 
dw 



or 



_ = _ Ve'"(l xe-««+^) -2iVe~* sin(a+ 0) 

<*J = - 7e-i,[( e ««+,)_ e -«a+,) )+ 2i s in(a + 0)] 
= -2iVe- i *[sin(a + <p) +sin(a+ 0)] 



and 



dw 



gj = 2 F[sin(a + ip) + sin(«+ £)] 
In general, for arbitrary r we find that 

we start with relation (10) 



To evaluate 



dz 



(36) 



(36') 



dz' 



CO p 

z' = ze lZ "- 



At the boundary surface 

z' = ze*-*°- l < where e = <p-8 and z = ae *<* +tr 



dz zV dz , 



(37') 



dp de A\ff ^ _ ^<fy 



dip 

di 



1+ de 



Then 

By equation (5) 
a 2 



dz; 

dz 



\de) 



i+— 
i+ de 



(37) 



f = s'+-> and at the boundary s' = ae*" He , or 
f = 2a cosh(^ + ii9) 



dz 



Then 



and 



= 2 sinh(^ + ^)e-^+ f «). 
df 



dz' 



= 4te~ 2 * (sinhV cos 2 0 + coshV sin 2 S) 
=4e~ 2 * (sinhV + sin 2 e) 



dr 



dz 

Then finally 



= 2e-*VsinhV + sin 2 0) 



(38) 



dw 




dw 




dz\ 


dz'l 


dr 




dz 




a?r 


drl 



F[sin(<* + <p) + sin (a + + g|)e*> 



(39) 



In this formula the circulation is given by equation 
(35). In general, for an arbitrary value of r (see 
equation (36')), the equation retains its form and is 
given by 



(40) 



(41) 



fo<« + '> + 4wI 1+ ifiV « 

■yj (sinAV + sin 2 0)(l + 
For the special case r = 0, we get 
Fsin(a + p)(l + j|Vo 

y(siDh 2 ^ + sin 2 9)(l+(^J) 

Equation (40) is a general result giving the velocity 
at any point of the surface of an arbitrary airfoil sec- 
tion, with arbitrary circulation for any angle of attack 
a. Equation (39) represents the important special 
case in which the circulation is specified by the Kutta 
condition. The various symbols are functions only of 
the coordinates {x, y) of the airfoil boundary and ex- 
pressions for them have already been given. In Tables 
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I and II are given numerical results for different air- 
foils, and explanation is there made of the methods of 
calculation and use of the formulas developed. 

We have immediately by equation (3) the value of 
the pressure p at any point of the surface in terms of 
the pressure at infinity as 

P 



5-'-(f)' 



Some theoretical pressure distribution curves are given 
at the end of this report and comparison is there made 
with experimental results. These comparisons, it will 
be seen, within a large range of angles of attack, are 
strikingly good. 18 

GENERAL WING-SECTION CHARACTERISTICS 

The remainder of this report will be devoted to a 
discussion of the parameters of the airfoil shape affect- 
ing aerodynamic properties with a view to determining 
airfoil shapes satisfying preassigned properties. This 
discussion will not only furnish an illuminating sequel 




Figure 16. 

to the foregoing analysis leading to a number of new 
results, but will also unify much of the existing theory 
of the airfoil. In the next section we shall obtain 
some expressions for the integrated characteristics of 
the airfoil. We start with the expressions for total 
lift and total moment, first developed by Blasius. 

Blasius' formulas. — Let C in Figure 16 represent a 
closed streamline contour in an irrotational fluid field. 
Blasius' formulas give expressions for the total force 
and moment experienced by C in terms of the complex 
velocity potential. They may be obtained in the fol- 
lowing simple manner. We have for the total forces 
in the x and y directions 

Px= - fpAs= - fpdy 
c c 



Py^fPv 
C 



fpdx 

c 



P x —iP v = ~ fp{dy+idx) 



» A paper devoted to more extensive applications to present-day airfoils is in 
progress. 



The pressure at any point is 

P=P° 
Then, 

P.-tP,-§/V(dy+»d*) 

* C 

ip r dwdv)y- 

~2{tedT Az 

where the bar denotes conjugate complex quantities. 
Since C is a streamline, vAy— v v dx — 0. Adding the 
quantity 

ipf (» v + Wx) fyxdy — v y dx) = 0 
c 



to the last equation, we get M . 



P,-iP,-%f(v.~ w y ) 2 (dx + idy) 



(42) 



The differential of the moment of the resultant 
force about the origin is, 
dMo = p(x dx+y dy) 

=i?. P. of p[xdx+y dy+i(ydx—xdy)] 
=B. P. of p z dz 
where "B. P. of" denotes the real part of the complex 
quantity. We have from the previous results 



Then 



dM 0 =-B. P. of 



2\dz) 



z dz 



and M 0 =-R. P. of | Jz dz 



(43) 



Let us now for completeness apply these formulas to 
the airfoil A in the f plane (fig. 14) to derive the Kutta- 
Joukowsky classical formula for the lift force. By 
equation (32) we have 

dz Ve 2*z + z* 
and by equation (25) 



df _«i_22_ 
dz 1 z 2 z 3 ■ 



Then 



dw_dw dz 
df "dz ' df 

- _ Ve*"- j£ - + (RWe-** - a, TV") - t + . . • 
2v z z 

» Cf. Blasius, H: Zs. f. Math. u. Phys. Bd. 58 S. 93 and Bd. 59 S. 43, 1910. 
Similarly, 

a less convenient relation to use than (42). 

Note that when the region about C is regular the value of the integral (42) remains 
unchanged by integrating about any other curve enclosing C. 



t 
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and 



/dw\ 2 _ 



z z 



where 



Ao=V 2 e 2t " 



A, = - 2 W + 2a, V 2 e 2i « - 



Then 

P*-iP„= 



=|(2«JL t ) 



Therefore 



P z = pW sin a 
P„ = pVT cos a 



and are the components of a force pVT which is per- 
pendicular to the direction of the stream at infinity. 
Thus the resultant lift force experienced by the airfoil is 

L = pVT (44) 

and writing for the circulation T the value given by 
equation (35) 

L*=4wBpV 2 sin (o+j3) (45) 

The moment of the resultant lift force about the 
origin f =0 is obtained as 

M.-S.P.of-f^)Vd f 
-«.P.of-|/(A-f + fH 



) 



X 



2 2 



=P. P. of-|2iri (coefficient of s _1 ) 
=#. P. of-|2«(^ 2 + 4,c,) 

or, M 0 is the imaginary part of tp(A 2 + AjcO- After 
putting 20 c x = me* and ai = 6 2 e 2<1, we get 

M> = 2irpV 2 b 2 sin 2 (a + 7) + pW m cos (a + h) (46) 

The results given by equations (44) and (46) have 
physical significance and are invariant to a transforma- 

» It-may be recalled that ci-- !? *Me'*ti v aud a, -a'+^'+ra. (See eq. (25').) 

ir o * 



tion of origin as may be readily verified by employing 
equations (26) and (32') and integrating around the d 
circle in the z x plane. It is indeed a remarkable fact 
that the total integrated characteristics, lift and loca- 
tion of lift, of the airfoil depend on so few parameters 
of the transformation as to be almost independent of 
the shape of the contour. The parameters B, p, a u 
and Ci involved in these relations will be discussed in a 
later paragraph. 

We shall obtain an interesting result 21 by taking 
moments about the point f = Ci instead of the origin. 
(M in fig. 17.) By equation (25) we have, 

and by equation (43) 

M M =B. P. 0 f-|/[^J( f - Cl )df 

=P.P.of-|/(^ + 4« + f + ...)x 

( f+ 5 + j + ...Xi-s + .> 

=R. P. of -wpA 2 



or 



M M = 2*b 2 pV 2 sin 2 (a + y) 




(47) 



x Axis 



Figure 17.— Moment arm from Monto the lift vector 

This result could have been obtained directly from 
equation (46) by noticing that pVT in the second term 
is the resultant lift force L and that Lm cos (a + 5) 
represents a moment which vanishes at M for all values 
of a. (In fig. 17 the complex coordinate of M is 
f = me K , the arm OH %s m cos (a + 5) .) The perpendic- 
ular h M from M onto the resultant lift vector is simply 
obtained from M M = Lh M , 
as 

, b 2 sin 2 (a + 7) 



z 2Bsia (a+ff) 



(48) 



The intersection of the resultant lift vector with the 
chord or axis of the airfoil locates a point which may 
be considered the center of pressure. The amount of 
travel of the center of pressure with change in angle 
of attack is an important characteristic of airfoils, 
especially for considerations of stability, and will be 
discussed in a later paragraph. 

!I First obtained by K. von Mises. (Reference 6.) The work of von Mises forms 
an elegant geometrical study of the airfoil. 
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The lift force has been found to be proportional to 
sin (ce + /9) or writing a+/8 = a I 



(49) 



where a x may be termed the absolute angle of attack. 
Similarly writing a+y=a 2 



M M =27rb 2 pV 2 wa.2a 2 



(50) 



With von Mises (reference 6, Pt. II) we shall denote 
the axes determined by passing lines through M at 
angles j3 and y to the x axis as the first and second axes 
of the airfoil, respectively. (Fig. 1#.) The directions 
of these axes alone are important and these are fixed 
with respect %o a given airfoil. Then the lift L is 
proportional to the sine of the angle of attack with 
respect to the first axis and the moment about M to 



If this moment is to be independent of a, the coeffi- 
cients of sin 2a and cos 2a must vanish. 
Then 

b 2 cos 2y =Rr cos (j3 + a) 

and 

6 2 sin 2y=Br sin(j3 + <r) 

Hence, 

h s 

r—^and <r=2y— /3 

Then if we move the reference point of the moment to 

b 2 

a point F whose radius vector from M is ^« 2t_? > * Ge 

moment existing at F is for all angles of attack con- 
stant, and given by 



M P =2xp6 2 F 2 sin 2(7- 0) 



(51) 



0 7 R " r Direction 

1 *■ 1 —n of sc 

, f T 



axis 



thruM 




°ff/ow 



Figure 18.— Illustrating the geometrical properties of an 

the sine of twice the angle of attack with respect to 
the second axis. 

From equation (47) we note that the moment at any 
point Q whose radius vector from M is re u , is given by 

M Q = 2rpb 2 V 2 sin 2 (a + y) - Lr cos (a + o-) 

Let us determL.,^ whether there exist particular 
values of r and <r for which M Q is independent of the 
angle of attack a. Writing for L its value given by 
equation (45), 

M Q - 2r P b 2 V 2 sin 2(a 4- 7) - IvpRrV 2 sin (a + 0) cos (a + a) 
And separating this trigonometrically 

M a *=2irpV 2 [{b 2 cos 27- Er cos (/3+<r)) sin 2a 
+ (ft 2 sin 2y—Rr sin (0+ a)) cos 2a 
-Rr sin (0-<r)] 



airfoil (axes and lift parabola of the E. A. F. 19 airfoil) 

It has thus been shown that with every airfoil pro- 
file there is associated a point F for which the moment 
is independent of the angle of attack. A change in 
lift force resulting from a change in angle of attack 
distributes itself so that its moment about F is zero. 

From equation (47) it may be noted that at zero lift 
(i. e., a = — j8) the airfoil is subject to a moment couple 
which is, in fact, equal to My. This moment is often 
termed "diving moment" or "moment for zero lift." 
If Mr is zero, the resultant lift force must pass through 
F for all angles of attack and we thus have the state- 
ment that the airfoil has a constant center of pressure, 
if and only if, the moment for zero lift is zero. 

The point F, denoted by von Mises as the focus of 
the airfoil, will be seen to have other interesting prop- 
erties. We note here that its construction is very 
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simple. It lies at a distance ^ from M on a line making 

angle 27-/3 with respect to the x axis. From Figure 
18 we see that the angle between this line and the first 
axis is bisected by the second axis. 

The arm h F from F onto the resultant lift vector L 
(h F is designated FT in Figure 18; note also that FT; 
being perpendicular to L, must be parallel to the direc- 
tion of flow; the line TV is drawn parallel to the first 
axis and therefore angle VTF=a+ 0) is obtained as 
, jM F _ -b 2 sin 2(0 -y) 
flF ~ L ~ 2jBsin( a +/3) 
b 2 

or setting ^~2R sm 2 — 

hF "~sm(a+0) (52) 
But h F is parallel to the direction of a, and the relation 
h= — h F sin (a+ 0) states then that the projection of 
h F onto the line through F perpendicular to the first 
axis is equal to the constant h (h is designated FV in 
the figure) for all angles of attack. In other words, 
the pedal points T determined by the intersection of 
h F and L for all positions of the lift vector L lie on a 
straight line. (The line is determined by T and V in 
fig. 18.) The parabola is the only curve having the 
property that pedal points of the perpendiculars 
dropped from its focus onto any tangent lie on a 
straight line, that line being the tangent at the vertex. 
This may be shown analytically by noting that the 
equation of L for a coordinate system having F as 
origin and FV as negative x axis is 

x sin ai + y cos a, = A F = ~ sin ^ +0) 

By differentiating with respect to a, = a + 0 and elim- 
inating on we get the equation of the curve which the 
lines L envelop as y 2 = 4h(x + h). From triangle FVS 
in Figure 18, it may be seen that the distance 
b 2 

MF= ^ is bisected at S by the line TV; for, since 

FV=h = 2ll sin2 (f-/3) and angle .rW= 2(0-y), then 
b 2 

SF= ^g- It has thus been shown that the resultant lift 

vectors envelop, in general, a parabola whose focus is 
at F and whose directrix is the first axis. The second 
axis and its perpendicular at M, it may be noted, are 
also tangents to the parabola being, by definition, the 

It 

resultant lift vectors for a= — y and a = ^ — y, respec- 
tively. 

If the constant h reduces to zero, the lift vectors 
reduce to a pencil of lines through F. Thus a constant 
center of pressure is given by h = 0 or sin 2(0— y) =0 
which is equivalent to stating that the first and second 
axes coincide. The lift parabola opens downward 
when the first axis is above the second axis (0^>y); it 
reduces to a pencil of lines when the two axes are 



coincident (0 = y) and opens upward when the second 
axis is above the first (0<^y). 

W. Miiller 22 introduced a third axis which has some 
interesting properties. Defining the complex coordi- 
nate f 0 as the centroid of the circulation by 

and using equations (25) and (32) one obtains 
f 0 -c 1 = x 0 + ii/ 0 

where 



b 2 



Xo ^ 2sux(a+0) lRsilla+ R shl {a + 2y)] 



1 



y ° = 2sm( a +0) [Bcosa -R cos (a + 27)] 



(53) 



The equation of the lift vector lines referred to the 
origin at M and x axis drawn through M is 

b 2 sin (a + y) 



x cos a — ysm a= 



2R sin (a+0) 



(•54) 



and it may be seen that the point (%, 2/ 0 ) satisfies this 
equation. The centroid of the circulation then lies on 
the lift vectors. By elimination of a from equation 
(53) one finds as the locus of (xq, j/o) 



b 



2xo[R cos jS--g-cos (0-2y)] + 2y„[R sin 0 
+|sin (0-2y)]=R 2 -£ 2 



(55) 



which is the equation of a line, the third axis, and 
proves to be a tangent to the lift parabola. ■ Geomet- 
rically, it is the perpendicular bisector of the line FF' 
joining the focus to the point of intersection of the 
first axis with the circle. (Fig. 18.) 

The conformal centroid of the contour. — It has 
already been seen that the point M has special inter- 
esting properties. The transformation from the air- 
foil to the circle having M as center was expressed in 
the normal form and permitted of a very small e(<p) 
curve. (See p. 14.) It was also shown that the 
moment with respect to M is simply proportional to 
the sine of twice the angle of attack with respect to 
the second axis. We may note, too, that in the pres- 
entation of this report the coordinate of M, f = Ci 

r 2tt 

= -=- S $e iv d-<p, is a function only of the first harmonic 
T 0 

of the \[>(<p) curve. 

We shall now obtain a significant property of M 
invariant with respect to the transformation from air- 
foil to circle. We start with the evaluation of the 
integral 

ff||* 

« Reference 7, p. 169. Also Zs. fur Ang. Math. u. Mech. Bd. 3 S. 117, 1923. 

Airfoils having the same first, second, and third axes are alike theoretically in 
total lift properties and also in travel of the center of pressure, i. e., they have the 
same lift parabola. 
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where A is the airfoil contour, ds the differential of 
dz 



arc along A, and 



df 



, as will be recalled, is the magni- 



fication factor of the transformation f =/(z) mapping 
airfoil into circle; i. e., each element ds of A when 

12 gives dS the differential of arc in the 



magnified by 



df 



plane of the circle, i. e., |dz|. Then we have, 



/f 

A 



dz 
df 



ds= /f (2) |dz| and by equation (25), 
c 

.^+,+5+3+ . .)idzi 

= 2 J (c, +Re** + 1 e-'* + #, «-«* + 



S" 2< 



.)i2d? 



= 2ir.Rc 1 

- Cl /ds= dy- 
er X 



ds 



Then 



/f 


dz 
df 


ds 


f 

i 


dz 
df 


ds 



(56) 



The point M of the airfoil is thus the conformal- cen- 
troid obtained by giving each element of the contour 
a weight equal to the magnification of that element, 
which results when the airfoil is transformed into a 
circle, the region at infinity being unaltered. It lies 
within any convex region enclosing the airfoil contour. 23 

ARBITRARY AIRFOILS AND THEIR RELATION TO 
SPECIAL TYPES 

The total lift and moment experienced by the air- 
foil have been seen to depend on but a few parameters 
of the airfoil shape. The resultant lift force is com- 
pletely determined for a particular angle of attack by 
only the radius R and the angle of zero lift 0. The 
moment about the origin depends, in addition, on the 
complex constants C\ and a y or, what is the same, on 
the position of the conformal centroid M and the focus 
F. The constants Ci and ai were also shown (see foot- 
note 20) to depend only on the first and second har- 
monics of the «(*>) curve. Before studying these 
parameters for the case of the arbitrary airfoil, it will 
be instructive to begin with special airfoils and treat 
these from the point of view of the "conformal angular 
distortion" [e(<p)] curve. 

Flow about the straight line or flat plate. — As a 
first approximation to the theory of actual airfoils, 
there is the one which considers the airfoil section to 
be a straight line. It has been seen that the line of 
length 4a is obtained by transforming a circle of radius 

a 2 

a, center at the origin, by f = zH — The region ex- 

z 



» Cf. P. Frank and K. Lowner, Math. Zs. Bd. 3, S. 78, 1919. Also reference 5, 
p. 146. 



ternal to the line 4a in the f plane maps uniquely into 
the region external to the circle |z| =a. A point Q Of 
the line corresponding to a point P at ae w is obtained 
by simply adding the vectors a(e u +e~ ie ) or completing 
the parallelogram OPQP'. 

For ^=0, we have from equation (6) . 
x = 2a cosh i> cos 8 = 2a cos 0 
y=2a sinh ^ sin 0=0 

Then" the parameters for this case are R = a, 0=0, 

ai=a 2 (i. e., b=a, 7=0), and M is at the origin O. 

Taking the Kutta assumption for detennining the 

circulation we have, 

the circulation, T = AnaV sin a 
the lift, L = 4xapF 2 sin a 

moment about M, M m =2to, 2 pV 2 sin 2a 

b 2 

position of F is at z r = c x + -g e"- 2 ^ = a 

Since 0 = 7, we know that the travel of the center of 
pressure vanishes and that the center of pressure is at 



(57) 




Figure 19. 



For at one-fourth the length of the line from the lead- 
ing edge. The complex flow potential for this case is 

. w(f) - - V[z^)e < "+~ e-'-l+l? log z(f) (58) 

where z ^)~\—-^(^ +a ^ s mverse °* e q ua *i° n 
(5). Since ^(¥>) = e(*») = 0 for this case, equation (39) 
giving the velocity at the surface reduces to 



sin I 



sin 



2 J 



for r = 45raP r sin a, 



and by equation (41) v- ^ (^J^*^ ) f° r r = 0- 

Flow about the elliptic cylinder. — If equation (5) 
is applied to a circle with center at the origin and 
radii"* ae*, the ellipse (fig. 19) 

3? y 2 
■ (2a cosh*) 2+ (2a sinh*) 2=1 . 

is obtained in the f plane and the region external to 
this ellipse is mapped uniquely into the region external 
to the circle. The same transformation also trans- 
forms this external region into the region internal to 
the inverse circle, radius ae~*. We note that a point 
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Q of the ellipse corresponding to P at ae^ w is 
obtained by simply completing the parallelogram 
OPQP' (fig. 19) where P' now terminates on the circle 
ae~*. The parameters are obtained as B=ae*, /3=0, 
Oj = (H, M is at the origin O. Then, assuming the rear 
stagnation point at the end of the major axis, 

Y=4:Ttae*V sax a 
£=4xpac*F s sin a 
M JI =2Ta 2 pV 2 sin 2a 

Since j8=y, the point F is the center of pressure for all 
angles of attack and is located at z F =ae~* from 0 or a 
distance ae* from the leading edge. The quantity 



EF 
EE'" 



ae* cosh f + sinh ^ I,, , . , A 
' 2a(e*+e-*) 4^sh^ i( 1 + tanh ^ 



represents the ratio of the distance of F from the 
leading edge to the major diameter of the ellipse. 

The complex flow potential is identical with that 
given by equation (58) for the flat plate, except that 
the quantity a 2 in the numerator of the second term is 
replaced by the constant aV*. Since ip(<p) = constant, 
«(«>) = 0 and equation (39) giving the velocity at each 
point of the surface for a stagnation point at end of 
major axis becomes 



p _ ^sin (y 4- a) + sin a]e* 
-/sinh V+sinV 

and for zero circulation by equation (41) 



sin(y +a)e* 



Vsinh V+sinV 



(59) 



(59') 



Circular arc sections. — It has been shown that 



the transformation £= 



= zH — applied to a circle with 
z 

center at 2=0 and radius a gives a straight line in the 
£ plane, and when applied to a circle with center 2=0 
and radius different from a gives an ellipse in the f 
plane. We now show that if it is used to transform a 
circle wi th cen ter at z=* is (s being a real number) and 
radius t/cP+s 2 , a circular arc results. The coordinates 
of the transform of the circle C in the f plane are given 
by equation (6) as 

x=2a cosh i> cos 0 
y=2a cosh ^ sin 6 

A relation between ^ and 0 can be readily obtained. 
In right triangle OMD (fig. 20), OM= s, angle OMD=d, 
and recalling that the product of segments of any 
chord through O is equal to a 2 , OD = % (OP-OP0 = 



=a sinh Then s sin 6=a sinh ^, and from 



the equation for y,y=2s sin 2 0. 
if in equation (6) we get 



Eliminating both 8 and 



the equation of a circle; but since y can have only 
positive values, we are limited to a circular arc. In 
fact, as the point P in Figure 20 moves from A' to A 
on the circle, the point Q traverses the arc Ai Ai and 
as P completes the circuit AA' the arc is traversed in 
the opposite direction. As in the previous cases, we 
note that the point Q corresponding to either P or to 
the inverse and. reflected, point P' is obtained by com- 
pleting the parallelogram OPQP'. We may also note 



B Curve ae*"'* 




FlOTJBE 20.— The circular arc airfoil 



that had the arc A\Ai been preassigned with the 
requirement of transforming it into the circle, the most 
convenient choice of origin of coordinates would be 
the midpoint of the line, length 4a, joining the end 
points. The curve B then resulting from using trans- 
formation (5) would be a circle in the z' plane, center 
at z'=>is, and the theory developed in the report could 
be directly applied to this continuous closed B curve. 




x Axis 



Figure 21.— Discontinuous B curve 

Had another axis and origin been chosen, e. g., as in 
Figure 21, the B curve resulting would have finite 
discontinuities at A and A', although the arc AiAi is 
still obtained by completing the parallelogram OPQP'. 

The parameters of the arc AiAi of chord length 4a, 
and maximum height 2s are then, R= Va a +s 2 , 



fi^taxr 1 -- The 
a 



focus F may be constructed by 
erecting a perpendicular to the chord at A' of length s 
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and projecting its extremity on MA'. The center M' 
of the arc also lies on this line. 

The infinite sheet having the circular arc as cross 
section contains as a special case the flat plate, and 
thus permits of a better approximation to the mean 
camber line of actual airfoils. The complex flow poten- 
tial and the formulas for the velocity at the surface 
for the circular arc are of the same form as those 
given in the next section for the Joukowsky airfoil, 
where also a simple geometric interpretation of the 
parameters e and f are given. 

Joukowsky airfoils. — If equation (5) is applied to 
a circle with center at a— s, s being a real number, and 
with radius B=a+s, a symmetrical Joukowsky air- 
foil (or strut form) is obtained. The general Joukow- 
sky airfoil is obtained when the transformation 
a 3 

f=sH — is applied to a circle C passing through the 

point 2= —a and containing z= a (near the circum- 
ference-usually), and whose center M is not limited to 
either the x or y axes, but may be on a line OM inclined 
to the axes. (Fig. 22.) The parametric equations of 
the shape are as before 

x=2a cosh $ cos 01 

(6) 

y=2a sinh ^ sin 0 

Geometrically a point Q of the airfoil is obtained by 
adding the vectors ae* +u and ae~*~ u or by completing 
the parallelogram OPQP' as before, but now P' lies on 
another circle B' defined as z=aer*~ i> , the inverse 
and reflected circle of B with respect to the circle of 
radius a at the origin (obtained by the transformation 
of reciprocal radii and subsequent reflection in the x 
axis). Thus OP-OP' — a % for all positions of P, and 
OP' is readily constructed. The center Mi of the 
circle B' may be located on the line AM by drawing 
OMi symmetrically to OM with respect to the y axis. 
Let the coordinate of M be z=is+de tf , where d, s, 
and 0 are real quantities. The circle of radius a, with 
center M a at z= is, is transformed into a circular arc 
through A, Ay which may be considered the mean 
camber line of the airfoil. At the tail the Joukowsky 
airfoil has a cusp and the upper and lower surfaces 
include a zero angle. The lift parameters are 

R~j¥+?+d, ^=tan-^. a 1 = a a =6V'T or b=a and 

7=0. Since 7=0, the second axis has the direction of 
the x axis. The focus F is determined by laying off 

the segment MF=-g on the line MA'. This quantity, 

it may be noted, is obtained easily by the following 
construction. In triangle MDC, MD=B, MC and 
MO are made equal to a, then CF drawn -parallel to 

a 2 

DC determines MF=>jg The lift parabola may be 

now determined uniquely since its directrix AM and 
focus F are known. 



It may be observed that if it is desired to transform 
a preassigned Joukowsky profile into a circle, there 
exists a choice of axis and origin for the airfoil such 
that the inverse of transformation (5) will map the 
airfoil directly into a circle. This axis is very approx- 
imately given by designating the tail as (—2a, 0) and 
the point midway between the leading edge and the 
center of curvature of the leading edge as (+ 2a, 0) the 
origin then bisecting the line joining these points. 




Floras 22.— The Joukowsky airfoil 



The complex potential flow function for the Jou- 
kowsky airfoil is 

w(f) - ~ F[0(r)e<«+^]+f£log0G-) (61) 
where 

By equation (39) we have for the velocity at the 
surface 

F[sin(« + <p) + sin (a + 0)] (l + ^ e*» 

2? = 1 .. ^ ■ „ 

y(sinhV+sin 2 0)(l+(^J) 

This formula was obtained by transforming the flow 
around C into that around B and then into that around 
A. Since we know that B is itself a circle for this 
case, we can simply use the latter two transformations 
alone. 
We get 

r y[sin(q+ y )+sin(a+g)]e* (62) 
Vsinh V+sin 2 0 

That these formulas are equivalent is immediately 
evident since the quantity 

is unity being the ratio of the magnification of each 
arc element of C to that of B. (See eq. (37).) 
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A very simple geometrical picture of the parameters 
« and ^, exists for the cases discussed. In Figure 
23 the value of t or at the point P is simply 




Figure 23.— Geometrical representation of < and * 
lor Joukowsky airfoils 

angle OPM, i. e., the angle subtended at P by the 
origin 0 and the center M. The angle of zero lift is 
the value of e for 8 = t; i e, t T m = (i = OTM. In 
particular, we may note that «=0 at S and /Si which 
are on the straight line OM. Consider the triangle 

OMP, where OP=ae*, MP=R=ue*°, jfp=f>, angle 

OPM=t; also, MOX=S, MOP = 0-5, OMP=x- 
(ip—h). Then by the law of cosines, we have 

= i + 2 P cos(#> - 5) + p 2 

or 



*->.=g log (1 + 2p cos(*>-«) +P 2 ) 



(63) 



and by the law of sines 



sin e 



p sin(y — d) 



or 



(l+2p cos(^-«) + p 2 )'/' 

/ \ L _i p sin(y-— 5) 
«(«>) = tan 't-^ — l s\ 



(64) 



We see that, as required, the expressions for the "radial 
distortion" and the "angular distortion" e(<p) 
are conjugate Fourier series and may be expressed as 
a single complex quantity 

-(_!)»-« 



(*-*o)-w=2- 



= -log[l + pe- << « , -»] 

It is evident also that the coefficient for n = l or the 
"first harmonic term" is simply pe ts and a translation 
by this quantity brings the circle C into coincidence 
with B as was pointed out on page 13. 



■, 'It 

The constant ^ 0 = s— fi>&<P is readily shown to be 



0 



invariant to the choice of origin 0, as long as 0 is 
within B. We have 

J- ^ 2 fhog (l+2p cos (,,-5)+pV*od¥> 



2ir 



0 



2 *0 2 




Figure 24.— The Joukowsky airfoil p=0.10, S-=45° 

Figure 24 shows the Joukowsky airfoil denned by 
p=0.10 and 5=45°, and Figure 25 shows the 
rp(<p), e(0), and e(*>) curves for this airfoil. 
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Figure 25.— The «(») and *(») curves for the airfoil in Figure 24 

Arbitrary sections.-* -In order to obtain the lift 
parameters of an arbitrary airfoil, a convenient choice 
of coordinate axes is first made as indicated for the 
Jbukowsky airfoil and as stated previously. (Page 7.) 
The curve resulting from the use of transformation (5) 
will yield an arbitrary curve ae* +i$ which will, in 
general, differ very little from a circle. The inverse 
and reflected curve ae~*~ u will also be almost circular. 
The transition from the curve ae* +it> to a circle is 
reached by obtaining the solution *•(*>) of equation 
(13). The method* of obtaining this solution as 
already given converges with extreme rapidity for 
nearly circular curves. 
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The geometrical picture is analogous to that given 
for the special cases. In Figure 26 it may be seen that 
a point Q on the airfoil (N. A. C. A. -M6) corre- 
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Figuee 28.— The N. A. C. A. -M6 airfoil 

sponding to P on the B curve (or P' on the B' curve) 
is obtained by constructing parallelogram OPQP'. 
The 4>(d) and «(0) curves are shown in Figure 27 for 
this airfoil. The complex velocity potential and the 
expression for velocity at the surface are given respec- 



The method used for arbitrary airfoils is readily 
applied to arbitrary thin arcs or to broken lines such 
as the sections of tail surfaces form approximately. In 
Figure 26 the part of the airfoil boundary above the x 
axis transforms by equation (5) into the two discon- 
tinuous arcs shown by full lines, while the lower 
boundary transforms into the arcs shown by dashed 
lines. If the upper boundary surface is alone. given 
(thin airfoil) we may obtain a closed curve ae* +u only 
by joining the end points by a chord of length 4a and 
choosing the origin at its midpoint. 25 The resulting 
curve has two double points for which the first deriva- 
tive is not uniquely denned and, in general, it may be 
seen that infinite velocities correspond to such points. 

At a point of the $(B) curve corresponding to a 
mathematically sharp corner, there exist two tangents, 



that is, the slope ^jjp is finitely discontinuous. 



The 
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Figure 27.— The and *(») curves for the N. A. C. A. -M6 airfoil 



tively by equations (33) and (39). The lift param- 
eters are 

R T 

R = ae* 0) p = e T aii (at0=ir),Mis atz = Ci= — f iK^e'^d? 

* 0 

and F is at z = c t + ^ where ai is given in equation (25')- 

The first and second axes for the N. A. C. A. -M6 
airfoil are found to coincide and this airfoil has then a 
constant center of pressure at F. Figures 28 (a) to 
28 (1) give the pressure distribution (along the x axis) 
for a series of angles of attack as calculated by this 
theory and as obtained by experiment. 2 * Table I 
contains the essential numerical data for this airfoil. 

« The experimental results are taken from test No. 323 of the N. A. C. A. variable- 
density wind tunnel. The angle of attack « substituted in equation (39) has been 
modified arbitrarily to take account of the effects of finite span, tunnel-wall inter- 
ference, and viscosity, by choosing it so that the theoretical lift is about 10 per cent 
more than the corresponding experimental value. The actual values of the lift 
coefficients are given in the figures. 



curve e(0) must have an infinite slope at such a point 
for according to a theorem in the theory of Fourier 
series, at a point of discontinuity of a F. S., the con- 
jugate F. S. is properly divergent. This manifests 
itself in the velocity-formula equation (39) in the fac- 
tor 1^1 + which is infinite at these sharp corners. 

For practical purposes, however, a rounding of the 
sharp edge, however small, considerably alters the slope 

djjS at this point. 

Ideal angle of attack. — A thin airfoil r represented 
by a line arc, has both a sharp leading edge and a 
sharp trailing edge. The Kutta assumption for fixing 
the circulation places a stagnation point at the tail for 
all angles of attack. At the leading edge, however, 

» Note that - -7M for this case. 
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Fiora is 28 a to e.— Theoretical and experimental pressure distribution lor the M6 airfoil at various angles of attack 
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Figures 28 f to J.— Theoretical and experimental pressure distribution for the M6 airfoil at various angles ol attack 
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the velocity is infinite at all angles of attack except 
one, namely, that angle for which the other stagnation 
point is at the leading edge. It is natural to expect 
that for this angle of attack in actual cases the Mo- 
tional losses are at or near a minimum and thus arises 
the concept of "ideal" angle of attack introduced by 
Theodorsen (reference 8) and which has also been 
designated "angle of best streamlining." The defini- 
tion for the ideal angle may be extended to thick 
airfoils, as that angle for which a stagnation point 
occurs directly at the foremost point of the mean 
* camber line. 

The lift at the leading edge vanishes and the change 
from velocity to pressure along the airfoil surface is 
usually more gradual than at any other angle of attack. 

Theoretical Experimental 
O Upper surface 
x Lower surface 
^ {Average flN. = 5 x10 s ) 

-4.0 - 



of this function, one can determine airfoil shapes of 
definite properties. The e(?>) function, which we have 
designated conformal angular distortion function, will 
be seen to determine not only the shape but also to 
give easily all the theoretical aerodynamic character- 
istics of the airfoil. 
An arbitrary e(<p) curve is chosen, single valued, of 

period 2t, of zero area, and such that — <» gj S 1. 

These limiting values of ^ are far beyond values 

yielding airfoil shapes. 57 The \^(<p) function, except for 
the constant <A>, is given by the conjugate of the 
Fourier expansion of e(<p) or, what is the same, by 
evaluating equation (14) as a definite integral. The 



Theoretical 




Experimental 
O Upper surface 
x Lower surface 
(Average R./V.=5x IQ S ) 



a = 18° 




ao 1 - 

Figtjres 28 k to 1.— Theoretical and experimental pressure distribution for the M6 airfoil at various angles ol attack 



The minimum profile drag of airfoils actually occurs 
very close to this angle. At the ideal angle, which we 
denote by a z , the factor [sin (a+p) + sin (a +&)] in 
equation (38) is zero not only for 0 = ir or« = « r =/3 but 
also for 0 — 0 or e = We get 

<*i+*tf— — («/+<r) or 

ai= -(j»+lll (65) 

CREATION OF FAMILIES OF WING SECTIONS 

The process of transforming a circle into an airfoil is 
inherently less difficult than the inverse process of 
transforming an airfoil into a circle. By a direct appli- 
cation of previous results we can derive a powerful and 
flexible method for the creation of general families of 
airfoils. Instead of assuming that the f(d) curve is 
preassigned (that is, instead of a given airfoil), we 
assume an arbitrary if/(<p) or t(<p) curve 28 as given. 
This is equivalent to assuming as known a boundary- 
value function along a circle and, by the proper choice 

« Subject to some general restrictions given in the next paragraph. 



constant ^ 0 is an important arbitrary 28 parameter 
which permits of changes in the shape and for a certain 
range of values may determine the sharpness of the 
trailing edge. 

We first obtain the variable 6 as 6 {<p) = <p — e (<p), so 
that the quantity ^ considered as a function of 9 is 
$(6) — $ [<p (6)]. The coordinates of the airfoil surface 
are then 

x=2a cosh 4> cos d 

(6) 

y=2a sum f sin 8. 

" For common airfoils, with a proper choice of origin, ||^| «0.30. 

a For common airfoils it is usually between 0.05 and 0.1S. The constant is 
not, however, completely arbitrary. We have seen that the condition given by 
equation (22) is sufficient to yield a contour free from double points in the 2' plane. 
We may also state the criterion that the inverse of equation (5) applied to this 
contour shall yield a contour in the f plane free from double points. Consider the 
function f(9) for 9 varying from 0 to » only. The negative of each value of *(9) in 
this range is considered associated with —9. i. e., x£9<2x. Designate the function 
thU3 formed from 9=0 to 2* by J (9)*. Then J(9)* represents a line arc in the f 
plane, i. e„ the upper surface of a contour. [See footnote 25.] Then for the entire 
contour to be free from double points it is necessary that the lower surface should not 
cross the upper, that is, the original f(9) curve for 9 varying from «• to 2r must not 
cross below f(9)\ 
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The velocity at the surface is 



V 



[sin (a+ <p) +sin (« + <3)] e+° 



(sinh V + sin 2 



(39') 



and is obtained by using equation (37') instead of (37) 
in deriving (39). The angle of zero lift /3 is given by 
<p(d)=d+i(e) for0 = x,i. e., <p (*■) =t+ 0. 
The following figures and examples will make the 
process clear. We may first note that the most natural 
method of specifying the e (<e) function is by a Fourier 
series expansion. In this sense then the elementary 
types of 6 (<p) functions are the individual terms of 
this expansion. 



29(h) to 29 (t). In particular, the second harmonic 
term may yield S shapes, and by a proper combina- 
tion of first and second harmonic terms, i. e., by a 
proper choice of the constants A u A 2 , S u and 5 2 in the 
relation 

e(<p) =Ai sin (<p— Si) + A 2 sin (2<p — 82) 

it is possible to fix the focus F of the lift parabola as 
the center of pressure for all angles of attack. 29 The 
equation 

e((e) =0.1 sin (<e-60°) + 0.05 cos 2<p 
represents such an airfoil and is shown in Figure 29 (u). 

The general process will yield infinite varieties of 
contours by superposition of sine functions; in fact, if 



e(ip)=0.l sin<p 

W : 


e{pJ-0.05 sinfp-45'J 
^ p^ a05 


e(pJ=-OJ sin 3p 
00=0-1 


tfoj-a/ sin(<p-l5°) 

ri — - — 


t/pJ-O./ sin(p-45') 
d/ 0 '0.08 

^ (i) ' 


e(p)-O.I sin(3p-45°) 


efpJ'O./ sin(p-30°) 
tpo-0.1 


efpj-a/ sin(p-45°) 

_>a- 


t(p)-O.I sin(3p-90°) 
— — 0o-Oi_^ 


e{p)-O.I sin(p-45°) 


t(ip)=O.I5 sin(p-45°) 




e(p)-0.075 sin(4pj 
--^ 00-0.075 


efpJ-OJ sin(<p-60°) 

0.-O./ ____ 


efpJ-O.i sin Bp 


— ^cs» " g:: --^ mSB ~ 

t(p)-0.075 sin(4p-45°) 


efpJ=OJ sinfp- 75') 

0O-O-I ■ 


t(p)-O.I sin(2p-45°) 
0.-O.I 

---^ W j _ - 


e(p)=0.075 sin(4p-90°) 
00-0.075 



€(p)-O.I sin(p-90°) 

00=0.1 



efpJ-OJ sin(2<p-90°) 

Figure 29.— Airfoils created by varying iM 




t(p)-O.I sin(p-60°)-.05 sin(2p-90°) 

0o=O.I 

(u) 



Consider first the effect of the first harmonic term 
c (<p) =Ai sin (<p-S0, &> = c 
In Figures 29(a) to 29(g) may be seen the shapes 
resulting by displacing 5i successively by intervals of 
15° and keeping the constants .At = 0.10 and ^0=0.10. 
The first harmonic term is of chief influence in deter- 
mining the airfoil shape. The case «(v) = 0.1 sin 
(p -45°)„is given detailed in Table II. (This airfoil 
is remarkably close to the commonly used Clark Y 
airfoil.) The entire calculations are characterized by 
their simplicity and, as may be noted, are completely 
free from the necessity of any graphical evaluations or 
constructions. 

The effect of the second and higher harmonics as 
well as the constant fa may be observed in Figures 



the process is thought of as a boundary-value problem 
of the circle, it is seen that it is sufficiently general to 
yield every closed curve for which Riemann's theorem 
applies. 



Langley Memorial Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., November 4> 193$. 

" This is accomplished as follows: We seek to determine the constants Au A j, Si, and 



tal; however, with but a few practice trials, solutions can be obtained at will. Addi- 
tion of higher harmonics will yield further shapes having the same center of pressure 
properties if f> is kept unchanged. 
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I. EVALUATION OF THE INTEGRAL. 



2x 



«(*') --f v f *M cot d v 



jr "q d<p 



sin 



dp 



(13) 
(13') 



The function if(<p) is of period 2x and is considered 
known. (Note that the variables <p and <p' are re- 
placed by 8 and 8', pi and p/, ^ and ', etc., in 
equation (21) and that the following formula is 
applicable for all these cases.) 

A 20-point method for evaluating equation (13) as 
a definite integral gives 

e Ws-J[«,¥ +aiWi- *-i) +<h(h- iM 



+ +a 9 (^,-^_ 9 )J,, 



a) 



where 



4i — value of 4(.<p) at ? = <p' + 



10 



^„ = value of 4(<p) at <ff = ^' + 



10 



(71 = 1,-1, 2, -2, 



. 9, -9). 



and the constants o„ are as follows: a o = Jo = 0.3142; 

Oi = 1.091; 02 = 0.494; 03 = 0.313; o 4 = 0.217; o 5 = 0.158; 
a 8 = 0.115; Or = 0.0884; a 8 = 0.0511; and 09 = 0.0251. 

This formula may be derived directly from the 
definition of the definite integral. The 20 intervals 1 

chosen are <p — ^ to <P+^fi' *° e *°' 

It is only necessary to note that by expanding 4(v) in 
a Taylor series around <p = <p' we get 

i S " m cot £^d^-2s 

<p'-s 2 L *v l<p=<p' 

where the interval <p' — s to <p' + s is small. And, in 



i f 4>M cot ^ d*> 



is very nearly 



- 4*a log 



sin 


<c 


— <P2 




2 


sin 




-<Pi 




2 



' Reference 2, p. 11, gives a 10-point method result. 
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where the range <P2 — <pi is small and 4a. is the average 
value of 4(<P) hi this range. The constants a„ for the 
20 divisions chosen above are actually 



smir- 



o„ = log 



271+1 



40 



sinir- 



2n-l 



40 



(ri= -9, 



+ 9) 



As an example of the calculation of e(8) we may refer 
to Table I and Figures 26 and 27 for the N. A. C. A. 
-M6 airfoil. From the 4/(8) curve (fig. 27) we obtain 

the 20 values of 4 and for 20 equal intervals of 8. 
For the airfoil (fig. 26) we get the following values: 



(Upper 
B surface) 




d? 

de 


(Lower 
9 surface) 




dj 

de 


0 (nose) 


0. 192 


0.000 


llx 
10 


0.049 


-0. 002 


X 

10 


. 185 


.027 


12x 
10 


.057 


.050 


2x 

10 


. 192 


.000 


13x 
10 


.071 


.030 


3x 
10 


. 189 


-. 030 


14x 
10 


.077 


.011 


4x 
10 


. 174 


-. 064 


15x 
10 


. 079 


. 000 


5x 
10 


.146 


-.095 


16x 
10 


. 082 


.016 


6x 
10 


. 110 


-. 114 


17x 
10 


. 090 


.039 


7x 
10 


.077 


-.086 


18x 
10 


. Ill 


. 091 


8x 
10 


. 052 


-.066 


19x 
10 


.150 


. 154 


9x 

10 


.041 


.025 


2x (nose) 


. 192 


. 000 


a- (tail) 


.055 


.000 









The value of e at the tail (i. e., the angle of zero Eft) 
is, for example, using formula I 

•=-;[fo*° 

+ 1.091 (.049 -.041) 
+ .494 (.057 -.052) 
+ .313 (.071 -.077) 
+ .217(.077-.110) 
+ .158(.079-.146) 
+ .115 (.082 -.174) 
+ .0884(.090-.189) 
+ .0511(.111-.192) 
+ .0251(.150- .185)] = .0105 
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3t 



The v alue of c f or 0 = > for example, is obtained by 
a cyclic rearrangement. Thus, 



".[lO 



(-.030) 



+ 1.091 (.174 -.192) 
+ .494(.146-.185) 
+ .313(.110-.192) 
+ .217(.077-.150) 
+ .158(.052-.lll) 
+ .115(.041-.090) 
+ .0884 (.055 -.082) 
+ .0511 (.049 -.079) 
+ .0251 (.057 - .077)] = .0347 

The 20 values obtained in this way form the ei(0) 
curve, which for all practical purposes for the airfoil 
considered, is actually identical with the final e(0) 
curve. 



II. NOTES ON THE TRANSFORMATION. 



r=/( 2 )=c 1 + z +| 1 +J+ 



(4') 



There exist a number of theorems giving general 
limiting values for the coefficients of the transforma- 
tion equation (4), which are interesting and to some 
extent useful. If {■=/(«) transforms the external 
region of the circle C of radius B in the z plane, into 
the external region of a contour A in the f plane in a 
one-to-one conforms! manner and the origin f=0 lies 
within the contour A (and/'(») = i) then the area S 
inclosed by A is given by the Faber-Bieberbach 
theorem as 2 



S = B?ir- 2 
71=1 



n\a„r 



Since all members of the above series term are positive, 
it is observed that the area of C is greater than that 
inclosed by any contour A in the f plane (or, at most, 
equal to the area inclosed by A if A is a circle). 
This theorem leads to the following results 

lail^B 2 (a) 
(4l£2B (b) 

Let us designate the circle of radius B about the 
conformal centroid M as center as Ci (i. e., the center 
is at f ~Ci; this circle has been called the "Grund- 
kreis" or "basic" circle by von Mises). Then since 

represents the distance of the focus F from M, uUe 

relation (a) states that the focus is always within &. 
In fact, a further extension shows that if r 0 is the radius 
of the largest circle that can be inclosed within A, then 

F is removed from & by at least • 

> For details of this and following statements see reference 5, p. 100 and p. 147, ana 
also reference 6, Part II. 



From relation (b) may be derived the statement that 
if any circle within A is concentrically doubled in radius 
it is contained entirely within a circle about M as 
center of radius 2B. Also, if we designate by c the 
largest diameter of A (this is usually the "chord" of 
the airfoil) then the following limits can be derived: 

B>\c 

These inequalities lead to interesting limits for the 
lift coefficient. Writing the lift coefficient as 

L 



%f>V*c 



where by equation (45) the lift force is given by 
Z = 4TBpF 2 sin (a+j3) 

we have 

8tB 



2Twa(.a+p)<,C L -- 



; sin(a+/3)^47Tsin(a+/3) (II) 



The flat plate is the only case where the lower 
limit is reached, while the upper limit is attained for 
the circular cylinder only. We may observe that a 
curved thin plate has a lift coefficient which exceeds 
2t sin (<*+ 0) by a very small amount. In general, the 
thickness has a much greater effect on the value of 
the lift coefficient than the camber. For common 
airfoils the lift coefficient is but slightly greater than 
the lower limit and is approximately 1.1X2tt sin 
(.a+fi). 

Another theorem, similar to the Faber-Bieberbach 
area theorem, states that if the equation f =/(z) trans- 
forms the internal region of a circle in the z plane into 
the internal region of a contour B in the f plane in a 
one-to-one conformal manner and/'(0) = 1 (the origins 
are within the contours) then the area of the circle is 
less than that contained by any contour B. This 
theorem, extended by Bieberbach, has been used in an 
attempt to solve the arbitrary airfoil. 3 The process 
used is one in which the area theorem is a criterion as 
to the direction in which the convergence proceeds. 
Although theoretically sound, the process is, when 
applied, extremely laborious and very slowly con- 
vergent. It can not be said to have yielded as yet 
really satisfactory results. 

III. LOCATION OF THE CENTER OF PRESSURE FOR AN 
ARBITRARY AIRFOIL 

It is~of some interest to know the exact location of 
the center of pressure on the x axis as a function of the 
angle of attack. In Figure 30, O is the origin, M the 
conformal centroid, L the line of action of the lift 
force for angle of attack a. Let us designate the 

' Muller, W., Zs. t. angew. Math. u. Mech. Bd. 5 S. 397, 1926. 
HShndorf, F., Zs. f. angew. Math. n. Mech. Bd. 6 8. 26S, 1926. 
Also reference 5, p. 185. 
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intersection of L with the x axis of the airfoil as the 
center of pressure P. 

In the right AONM we have, 

0M= c^me^^Ai+iBt 

ON=m cos 5=Ai 

MN= m sin 5 = Bi 

and in right AJKM, £M-^-> 



sin a sin a 
hit 



sm a 



Then 

and NP =*KN tan a = h M sec a— Bi tan a 

By equation (48) 

a -Mjt 6 2 sin 2 (a+y) 
lhM ~ L ~2R sin(«+0) 

Then the distance from the origin to the center of 
pressure P is 



OP = 0N+ NP =Ai-Bi tan a 
W sin 2 (a+y) 
2#cos a sin (a+P) 



L iff vector 



(III) 




x Axis 



Figure 30— Center of pressure location on the x axis 

EXPLANATION OF THE TABLES 

Table I gives the essential data for the transforma- 
tion of the N. A. C. A. -M6 airfoil (shownin fig. 26) 
into a circle, and yields readily the complete theoretical 
aerodynamical characteristics. Columns (1) and (2) 
define the airfoil surface in per cent chord; (3) and (4) 
are the coordinates after choosing a convenient origin 
(p. 7); (5) and (6) are obtained from equations (7) 
and (8) of the report; (9) is the evaluation of equation 
(13) (see Appendix); (10) and (11) are the slopes, ob- 



tained graphically, of the f against 0, and £ against 
0 curves, respectively; (12) is given by 



^(smhV+sm^H-^) 3 ) 
1 2x 

where ^o=7r S 4> (<p) d*> and may be obtained graphi- 

cally or numerically; column (13) gives p=0+e. The 
velocity v, for any angle of attack, is by equation (39) 

v=Vk [sin (a + ¥>)+sin (a+/3)] 
and the pressure is given by equation (3). The angle 
of zero lift 0 is the value of e at the tail; i. e., the value 
of e for 6 = w. 

Table II gives numerical data for the inverse process 
to that given in Table I; viz, the transformation of a 
circle into an airfoil. (See fig. 29.) The function 
e(<p)=0.1 sin (<p— 45°) and constant ^ 0 =0.10 are 
chosen for this case. Then ^(<p) =0.1 cos (<p— 45°) 
+ 0.10. It may be observed that columns (11) and 

(12) giving the coordinates of the airfoil surface are 
obtained from equations (6) of the report. Column 

(13) is given by 

and the velocity at the surface is by equation (39') 
v=Vk [sin (a+*>)+sin («+£)] 
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TABLE I 
N. A. C. A. — M6 

UPPER SURFACE 
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HHER SU3 


*FACE 








0 


0 


2.037 


0.000 


0.000 


0.0373 


6.283 


a 192 


-0.0487 


0.000 


0.088 


6.280 


-2 


37 


1.25 


1.76 


1.986 


-.071 


.0125 


.0297 


6.078 


.172 


-.0781 


.133 


.120 


4.615 


-16 


21 


2.50 


2.20 


1.936 


-.089 


.0844 


.0234 


5.989 


.152 


-.0860 


.160 


.080 


3.525 


-21 


43 


5.0 


2.73 


1.838 


-. 110 


.173 


.0176 


5.855 


.132 


-.0882 


.133 


-.010 


2.510 


-29 


35 


7.5 


3.03 


1.734 


-.122 


.259 


.0144 


5.749 


.120 


-.0850 


.109 


-.045 


2.025 


-35 


28 


10 


3.24 


1.633 


-.131 


. .342 


.0128 


5.689 


.112 


-.0811 


.080 


-.057 


1.764 


-40 


24 


15 


3.47 


1.431 


-.140 


.494 


.0099 


5.505 


.099 


-.0723 


.069 


-.087 


1.466 


-48 


44 


20 


a 62 


1.229 


-.146 


.626 


.0085 


8.371 


.092 


-.0637 


.067 


-.067 


1.307 


-55 


64 


30 


3.79 


.826 


-.153 


.831 


.0070 


5.136 


.084 


-.0616 


.025 


-.052 


1.166 


-68 


39 


40 


3.90 


.421 


-.158 


.966 


.0066 


4.924 


.081 


-.0421 


.008 


^.036 


1.098 


-80 


17 


50 


3.94 


.017 


-. 159 


1.000 


.0063 


4.712 


.079 


-.0360 


.000 


-.029 


1.081 


-91 


59 


60 


3.82 


-.387 


-.164 


.963 


.0062 


4.818 


.0788 


-.0310 


.010 


-.013 


1.120 


-102 


S3 


70 


3.48 


-.791 


-.141 


.845 


.0058 


4.307 


.076 


-.0300 


.019 


.000 


1.211 


-114 


65 


80 


2.83 


-1. 195 


-.114 


.645 


.0050 


4.074 


.071 


-.0298 


.036 


-.011 


1.370 


-128 


14 


90 


1.77 


-1.599 


-.072 


.363 


.0035 


3.788 


.069 


-.0235 


.044 


-.040 


1.769 


—144 


16 


98 


1.08 


-1.801 


-.044 


.191 


.0025 


3. 594 


.060 


-.0140 


.020 


-.067 


2.368 


-154 


60 


100 


.26 


-2.003 


.000 


.000 


.0030 


3.142. 


.056 


.0106 


.000 


-.027 


19.83 


-179 


24 



TABLE II 

•(iri-0.1 sin (*>-45°) *„=»0.10 -0.0657=3° 47' 

UPPER ■ SURFACE 



*> 




8 




d< 
d* 


d* 

"Si 












• V 








c 






* 


cost I 


sini <// 


COS 0 


sin* 


X 


* 


Degrees 


Radians 




Radians 


Degrees 










2 


2 




0 


h* 

0.0000 


-0.0707 


0.0707 


O 1 

43 


a 1707 


0.0707 


0.0707 


1.0146 


a 1718 


0.9976 


0.0706 


1.0121 


0.0121 


6.3941 


5 


.0873 


-.0643 


.1516 


8 41 


.1766 


.0766 


.0643 


1.0156 


.1776 


.9886 


.1510 


1.0039 


.0268 


5. 1215 


10 


.1746 


-.0574 


.2319 


13 17 


.1819 


.0819 


- .0674 


1.0186 


.1826 


.9733 


.2298 


.9895 


.0420 


4.0960 


15 


.2818 


-.0500 


.3118 


17 62 


.1866 


.0866 


.0500 


1.0178 


.1877 


.9618 


.3068 


.9685 


.0876 


3.3602 


20 


.3491 


-.0423 


.3914 


22 26 


.1906 


.0908 


.0423 


1.0182 


.1918 


.9243 


.3816 


.9411 


.0732 


2.8421 


25 


.4363 


-.0342 


.4706 


28 87 


.1940 


.0940 


.0342 


1.0189 


.1962 


.8914 


.4632 


.9082 


.0888 


2.4704 


30 


.5236 


-.0259 


.5495 


31 29 


.1966 


.0966 


.0259 


1.0194 


.1979 


.8528 


.5223 


.8893 


.1034 


2.1892 


35 


.6109 


-.0174 


.6283 


36 0 


.1985 


.0985 


.0174 


1. 0198 


.1998 


.8090 


.5878 


,8250 


.1174 


1.9748 


45 


.7854 


.0000 


.7854 


45 0 


.2000 


.1000 


.0000 


1.0201 


.2013 


.7071 


.7071 


.7213 


.1423 


1.6889 


55 


.9599 


.0174 


.9425 


54 0 


.1985 


.0935 


-. 0174 


1. 0198 


.1998 


.5878 


.8090 


.5994 


.1616 


1.4709 


70 


1.2217 


.0423 


1.1794 


67 38 


.1906 


.0906 


-.0423 


1.0182 


.1918 


.3813 


.9244 


.3882 


.1773 


1.2859 


80 


1.3963 


.0674 


1.3389 


76 43 


.1819 


.0819 


-.0574 


1.0166 


.1826 


.2298' 


.9733 


.2336 


.1777 


1.2133 


90 


1.5708 


.0707 


1.6001 


86 57 


.1707 


.0707 


-.0707 


1.0146 


.1715 


.0706 


.9975 


.0716 


.1711 


1. 1717 


100 


L 7453 


.0819 


1.6634 


95 18 


.1674 


.0574 


-.0819 


1.0124 


.1581 


-.0924 


.9957 


-.0935 


.1874 


1.1586 


110 


1.9199 


.0908 


1.8293 


104 49 


.1423 


.0423 


-.0906 


1. 0101 


.1428 


-. 2557 


.9668 


-.2583 


.1381 


1. 1756 


125 


2.1817 


.0986 


2.0832 


119 22 


.1174 


.0174 


-.0985 


1.0069 


.1177 


-.4904 


.8715 


-. 4938 


.1028 


1.2727 


135 


2.3562 


.1000 


2.2582 


129 16 


.1000 


.0000 


-.1000 


1.0080 


.1002 


-.6329 


.7742 


-.6361 


.0776 


1.4088 


180 


2.6180 


.0966 


2.8214 


144 28 


.0741 


-.0259 


-.0966 


1.0028 


.0742 


-.8138 


.6812 


-.8161 


.0431 


1.8306 


180 


2.7925 


.0906 


2.7019 


184 48 


.0577 


-.0423 


-.0906 


1. 0017 


.0677 


-.9048 


.4258 


-.9063 


.0246 


2.4584 


170 


2.9671 


.0819 


2.8882 


186 19 


.0426 


-.0574 


-.0819 


1.0009 


.0426 


-.9873 


.2538 


-.9682 


.0108 


4.0496 


180 


3. 1416 


.0707 


3.0709 


175 57 


.0293 


-.0707 


-.0707 


1.0004 


.0293 


-.9975 


.0708 


-.9979 


.0021 


13.4411 



LOWER SURFACE 



0 


0.0000 


-0.0707 


0.0707 


4 3 


0.1707 


0.0707 


.0707 


1. 0146 


0. 1715 


a 9975 


0.0706 


1.0121 


a 0121 


8.3941 


-5 


-.0873 


-.0766 


-.0107 


-0 37 


.1643 


.0643 


.0766 


1.0135 


.1650 


.9999 


-.0108 


1.0134 


-.0018 


7.1236 


-10 


-.1748 


-.0819 


-.0926 


-5 18 


.1574 


.0574 


.0819 


1.0124 


.1581 


.9957 


-.0924 


1.0080 


-.0146 


6.3827 


-15 


-.2618 


-.0868 


-. 1782 


-10 2 


.1500 


.0500 


.0896 


1.0113 


.1606 


.9847 


-.1742 


.9958 


-.0283 


5.0338 


-20 


-.3491 


-.0908 


-.2585 


-14 49 


.1423 


.0423 


.0906 


1. 0101 


.1428 


.9668 


-.2557 


.9766 


-.0368 


3.9225 


-25 


-.4363 


-.0940 


-.3423 


-19 37 


.1342 


.0342 


.0940 


1.0090 


.1346 


.9420 


-.3357 


.9605 


-.0462 


3.1489 
2.6077 


-30 


-.8236 


-.0968 


-.4270 


-24 28 


.1259 


.0259 


.0986 


1.0079 


.1282 


.9102 


-.4142 


.9174 


-.0823 


-35 


-.6109 


-.0985 


-.5124 


-29 21 


.1174 


.0174 


.0988 


1.0069 


.1177 


.8716 


-.4901 


.8776 


-.0877 


2.2203 


-45 


-.7884 


-.1000 


-.6854 


-39 16 


.1000 


.0000 


.1000 


1.0050 


.1002 


.7742 


-.6329 


.7781 


-.0634 


1.7161 


-55 


-.9899 


-.0985 


-.8814 


-49 21 


.0826 


-.0174 


.0985 


1.0034 


.0827 


.6514 


-.7687 


.6636 


-.0627 


1.4168 


-70 


-1.2217 


-.0908 


-1. 1311 


-64 48 


.0577 


-.0423 


.0906 


1.0017 


.0577 


.4258 


-.9048 


.4265 


-.0522 


1.1660 


-80 


-1.3963 


-.0819 


-1.3144 


-76 19 


. .0426 


-. 0574 


.0819 


1.0009 


.0428 


.2535 


-.9673 


.2637 


-..0412 


1.0763 


-90 


-1. 5708 


-.0707 


-1.5001 


-88 57 


.0293 


-.0707 


.0707 


1.0004 


.0293 


.0706 


-.9975 


.0708 


-.0292 


1.0322 


-100 


-1.7483 


-.0574 


-1.6879 


-98 43 


.0181 


-.0819 


.0574 


1.0002 


.0181 


-. 1170 


-.9931 


-. 1170 


-.0180 


1.0270 


-110 


-1.9199 


-.0423 


-1.8776 


-107 3 


.0094 


-.0906 


.0423 


i.oooo 


.0094 


-.3021 


-.9533 


-.3021 


-.0090 


1.0616 


-125 


-2.1817 


-.0174 


-2.1643 


-124 1 


.0015 


-.0985 


.0174 


1.0000 


.0015 


-.5594 


-.8289 


-.6594 


-.0012 


1.2136 


-135 


-2.3562 


.0000 


-2.3562 


-135 0 


.0000 


-.1000 


.0000 


1.0000 


.0000 


-.7071 


-.7071 


-.7071 


.0000 


1.4209 


-180 


-2.6180 


.0259 


-2.6439 


-161 29 


.0034 


-.0968 


-.0259 


1.0000 


.0034 


-.8787 


-.4774 


-.8787 


-.0016 


2.1108 


-160 


-2.7925 


.0423 


-2.8348 


-162 25 


.0094 


-.0906 


-.0433 


1.0000 


.0094 


-.9533 


-.3021 


-.9833 


-.0028 


3.3801 


-170 


-2.9671 


.0574 


-3.0245 


-173 18 


.0181 


-.0819 


-. 0574 


1.0002 


.0181 


-.9932 


-. 1167 


-.9934 


-.0021 


8.6641 


-ISO 


-3.1418 


.0707 


-3.2123 


-184 3 


.0293 


-.0707 


-.0707 


1.0004 


.0293 


-.9975 


+.0708 


-.9979 


+.0021 


13.4411 




■ ^parallel 



. ' ■ 9«m: 1 ' . Positive 



. - , symbol resignation ^ol 



gym- Positive... 



Designa- 

H tionv"* 

.VR>T;'r S57/ 

> • 

Roll---'- 



linear 

: (ciompo-'i 
-neat, along 
axis)- 



Angular 




of' control surface (relative to" neutral 



D Uiameier . , ,. 

p Geometric pitch — . : . • 
' * ' p/D Pitch ratio - < < - 



(Indicate surface by proper, subscript.) 



^jl^Kwer, absolute coef|«ient;6 P =^^^^ : -■; 

G, >: Sp^rP'6wer coefficients^ 

^:*|^^|^ieja6| 




v,ri - ... ,,. } 



